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1  Program  Overview 


LRIR13RI08COR 


Title:  Investigating  Quantum  Data  Encrypted  Modulation  States 
Reporting  Period:  FY14 

Laboratory  Task  Manager:  Dr.  David  H.  Hughes 

Commercial  Phone:315-330-4122  DSN:587-4122 

E-Mail  Address:  david.hughes.16@us.af.mil 
AFOSR  Program  Manager:  Dr.  Arje  Nachman 

Research  Objectives:  Investigate  and  measure  modulated  optical  signals,  both  quantum  and 
classical,  for  Information  confidentiality,  integrity,  and  availability 

1.1  Introduction 

AFOSR  LRIR13RI08COR  supports  four  related  efforts,  all  primarily  concerned  with  aspects  of  optical 
communications.  However,  some  of  our  efforts  also  include  the  RF  spectrum  in  that  we  are  working 
toward  a  combined  optical/RF  link  operating  between  our  site  in  Rome  NY  and  Stockbridge  NY.  We 
leverage  AFRL  6.2  funds  as  well.  AFOSR  FY  2014  funds  were  largely  devoted  to  salary  for  government 
and  onsight  contractors.  Some  funds  have  been  used  for  equipment  purchases.  A  breakdown  of  the 
spending  is  given  in  Table  1  below. 

1.2  Technical  Work 

Section  2  reports  on  basic  research  into  noise  corrupting  phase  estimation  in  coherent  states  of 
light  in  the  1.5  micron  band.  Over  and  above  the  quantum  noise  of  light,  with  which  recent 
physical  cryptographic  systems  utilize  to  impart  randomness  to  a  ciphertext,  additional  noise 
sources  can  have  an  even  more  debilitating  effect  on  attempts  to  decrypt  ciphertext.  These 
latter  sources  of  noise  include  detection  system  thermal  noise  as  well  as  source  laser  phase 
drift  noise. 

Section  3  reports  on  a  basic  comparative  study  between  long  wavelength  infrared  (LWIR)  versus 
short  wavelength  infrared  (SWIR)  propagation  in  fair  and  inclement  weather.  Currently  the 
communication  bands  under  study  are  1.5  and  10  micron  bands.  Despite  the  huge 
infrastructure  devoted  to  the  1.5  micron  band,  work  in  the  10  micron  band  by  several  workers, 
including  AFRL/RITE  engineers,  indicates  that  propagation  in  the  10  micron  band  is  more  robust 
in  free  space  optics.  That  translates  to  greater  information  availability. 
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Section  4  reports  on  the  combined  optical/RF  link  work  under  way  between  Stockbridge  and 
Rome.  Currently,  the  optical  band  is  primarily  1.5  micron  band  carrier  frequencies;  turbulence 
mitigation  is  accomplished  by  curvature  adaptive  optics.  Planned  is  deployment  of  the 
LWIR/SWIR  comparison  experiment  system  over  the  longer  range  provided  by  the 
Rome/Stockbridge  link  as  well  as  a  coherent  state  quantum  data  encryption  system. 

Section  5  characterizes  the  propagation  of  entangled  photon  pairs  through  a  hyperspectral 
filter  device  originally  designed  for  multi-access  laser  communications  between  a  hub  and 
multiple  spokes.  We  asked  the  question:  Is  polarization  entanglement  preserved  in  transit 
through  the  device?  If  not,  can  the  device  be  modified  to  preserve  the  entanglement? 


Table  1.  Resource  Dispersal  FY14 


Government  Salaries 

$107,000 

Onsite  Contractor  Support 

$99,640 

Equipment 

$73,360 

Total 

$280,000 
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2  Coherent  State  Phase  Estimation 

(Mr.  Paul  Cook) 

2.1  Introduction  and  Background:  Expected  Results  via  Theory  and  Simulation 

The  estimation  of  phase  will  be  carried  out  by  heterodyne  and  homodyne  methods.  The  signal  will  be 


Figure  1.  An  Optical  90°  Hybrid  with  differential  detectors. 

mixed  with  a  local  oscillator  (LO)  using  a  90  degree  optical  hybrid  as  shown  in  Figure  1.  The  resulting  I 
and  Q.  signals  are  detected  with  differential  detectors  which  subtract  out  the  DC  component  and  any 
systematic  laser  amplitude  noise. 

We  can  write  the  signal  and  LO  as 


4sexp  (i(tcos  +  <ps) 

(2.1) 

Aloexp(i(ta)lo  +  cpl0 ) 

(2.2) 

where  AsanAAl0,  cos  and  &)to,  cps  and  (pio  are  the  detected  photon  amplitude  (a)  in  square  root 
photons  per  symbol  (i.e.,  in  a  digital  communication  system),  frequency,  and  phase  of  the  two  signals 
and  where  t  is  time.  The  relationships  between  the  amplitudes  and  the  powers  are 

As  =  J2PS  Al0  =  (2-3) 

The  signal  amplitude  inputs  and  outputs  for  the  second  two  splitters  on  the  right  in  Figure  1  are  shown 
n 

in  Figure  2  with  the  —  added  to  the  LO  signal  in  the  top  branch  (Q  branch). 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 

3 


“P  eXpO(tGX;  +  *S)) 

V2 


Alo 


e^i^to,0  +  <frlo  + 


Ai0  (f  * 

—  expl  I  tm,0  +  <h„  +  — 


Ah 


“p({‘0,I"  +  +1"  +  f)) 


2")) +  Tap^Uas 1  ♦»)) 

"Tesp({tos  +  *s+l)) 


Q 


-j='Mius>s+t$ 


e^>(i(ta|0  +  +lo)) 


•  ejqj^la^  1  +lu  1  Ajj  +  ~  exp(i(tG>s  +  +,)) 


+  *1.))  +  Y  e^i^tos  +  *s  +-  A  Jj 


Figure  2  The  final  signals  applied  to  the  photodetectors. 

Taking  the  real  part  of  each  of  the  four  terms  on  the  right,  squaring,  and  subtracting  the  top  two  from 
each  other  and  subtracting  the  bottom  two  from  each  other  (because  of  the  differential  detection),  we 
arrive  at  the  detected  signals  for  I  and  Q  by  extracting  the  frequency  difference  terms 


1 

I  =  -■  As  ■  Alo  ■  cos(A cp  +  A<x>t ) 

Q  =  ±-As-  Ai0  1  sin(A<p  +  A  cot)  (2.4a, b) 

/  =  sjPs  ■  Pi0  ■  cos(A <p  +  Acot) 

Q  —  yJ~Ps  ■  P/o  ■  sin(A<p  +  A  cot)  (2.5a,b) 


where  Ps  and  Plo  are  the  optical  power  of  the  signal  and  LO  respectively,  A  cp  is  the  relative  phase  of  the 
applied  signal  and  LO,  and  Ao)  is  the  frequency  difference  between  the  signal  and  the  LO. 


A  simulation  was  performed  in  MathCad  starting  with  the  four  equations  on  the  right  of  Figure  2.  In  each 
of  the  four  time  series  the  square  of  the  real  part  of  each  time  sample  was  taken  to  obtain  the  RF 
current  amplitude  (or,  equivalently,  modulated  detected  optical  power)  for  that  point  in  time  in  terms  of 
detected  photons  per  sample  point.  Each  point  was  then  "randomized"  by  assigning  it  a  value  by  the 
MathCad  Poisson  generator  to  generate  an  ensemble  for  each  of  the  four  detector  signals.  Ps  was  set  to 
25  (the  amplitude  then  is  V2  *  25)  detected  photoelectrons  and  Plo  was  set  to  1  x  10A6  (the  amplitude 
is  V2  *  106)  detected  photoelectrons  and  the  phase  angle  was  set  to  45  degrees.  Each  pair  of  resulting 
signals  for  the  I  and  Q.  branches  were  subtracted  to  simulate  the  current  subtraction  in  the  differential 
detection  diode  circuit.  A  Fourier  transform  was  performed  on  each  of  the  resulting  I  and  Q  signals  and 
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the  value  at  the  difference  frequency  for  each  was  recorded.  This  was  repeated  5000  times.  Each  of  the 
resulting  [I ,  Q]  values  is  plotted  below  in  Figure  3. 


v 


a 


u 


a 


Figure  3.  Heterodyne  simulation  of  I  (vg)  and  Q  (ug)  output  of  a  90  degree  optical  hybrid. 

Since  the  90  degree  hybrid  splits  the  signal  into  two,  the  I  and  Q.  branches  each  receive  25/2  signal 
photoelectrons  and  10A6/2  LO  photoelectrons  combined.  The  expected  signal  in  the  I  branch  when  the 
input  signal  is  at  a  phase  angle  of  45  degrees  and  is 

I  =  Vps  ■  Piocos(45) 

/  =^V25-10A6 

V2 

/  =  3535.5  (2.6) 

And  likewise  for  the  Q  branch 

Q  =  y[Ps  '  ho  sin  (45) 

Q  =  -j=V25  ■  10A6 

V2 

Q  =  3535.5  (2.7) 

for  a  total  vector  length  of 

R  =  V3535.52  +  3535. 52  =  5000  (2.8) 
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The  The  I  and  Q  branch  each  receive  ~  —  and  so  have  noise  each.  Added  in  quadrature  the  radial 

(amplitude)  noise  is  VZO  =  VlOA6  =  1000.  These  calculations  agree  well  with  the  simulation  within 
the  limits  of  the  noise  (Amplitude  =  5003.0  photons  /measurement,  amplitude  noise  =  1000.1  rms 
photons  /measurement) .  The  signal  to  noise  ratio  for  a  heterodyne  detected  signal  of  25  detected 
photons  using  a  90  degree  optical  hybrid  is  V25  =  5  which  also  agrees  well  (5003.0/1000.1  =  5.002)  with 
the  simulation  results. 

The  results  of  simulation  of  magnitude  and  phase  noise  of  the  system  as  described  was  performed  and  is 
shown  in  Figure  4.  This  simulation  was  carried  out  in  MathCad.  The  phase  noise  for  the  all  optical  full 
quadrature  detection  is 

( Pnoise  ~  ^  (2-9) 

For  n  =  25  we  should  have  a  phase  noise  of  0.2  Radians.  The  simulation  is  in  good  agreement  yielding 
0.204  R. 


Figure  4.  Simulation  of  phase  noise  in  the  heterodyne  90  hybrid  system.  The  theoretical  plot  for  (pnoise 


where  ny  is  the  number  of  photons  per  measurement  period. 
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Figure  4  plots  theory  and  simulation  of  phase  noise  vs  detected  signal  photons.  The  simulation  is 
noticeably  noisier  than  the  theory  for  very  low  photon  number  (<  10)  because  at  these  low  levels  the 
signal  amplitude  is  close  to  the  IQ  origin  and  outliers  of  the  simulated  measurements  are  in  different 
quadrants  than  that  of  the  centroid  of  the  signal  and  in  some  cases  there  are  simulated  data  points  in  all 
four  quadrants  as  shown  in  figure  5.  These  outliers  can  then  have  deviations  on  the  order  of  n/2 
whereas  the  theory  is  derived  from  a  small  angle  approximation  that  will  not  be  accurate  for  angles 
greater  than  0.1  radians. 
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Figure  5.  Here  n  =  2  (cp  =  45  degrees)  and  the  simulation  results  are  in  all  four  quadrants  producing  increased 

phase  noise  over  the  simplified  theoretical  model. 


2.2  Initial  measurements 

Measurements  have  been  made  with  a  simplified  system  shown  in  Figure  6  in  order  to  verify  the  basic 
heterodyne  system  and  scaling  factors  such  as  fiber  optic  connector  loss,  50/50  splitter  loss,  detector 
quantum  efficiency,  amplifier  gain,  and  noise  sources  over  and  above  photon  shot  noise.  The  signal  and 
Local  Oscillator  are  two  independent  lasers  attenuated  to  about  15  milliwatts  and  0.1  to  10  nanowatts. 
A  polarization  rotator  is  adjusted  to  align  the  linear  and  ellipsoidal  components  of  the  polarizations  of 
the  two  lasers  to  maximize  interference  in  the  50/50  splitter  (180  degree  hybrid).  The  heterodyne  signal 
is  differentially  detected  and  fed  to  a  low  noise  amplifier  before  being  digitized  by  a  Tektronix 
DP073304D  digital  Oscilloscope  at  100  Gigasamples/s.  For  these  tests  the  frequency  difference  between 
the  signal  and  LO  lasers  was  2.5  GHz  and  the  scope  lowpass  filters  were  set  at  4  GHz. 
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Oscilloscope 


Figure  6  Simplified  system  for  checkout 
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Figure  7  Shows  the  theoretical  and  simulated  phase  noise  of  the  180°  hybrid  with  the  90°  results 
included  for  comparison.  The  180°  hybrid  has  lower  noise  because  it  does  not  split  the  light  signal 
between  I  and  Q  branches  as  in  the  90°  hybrid.  With  the  180°  hybrid  the  l/Q  processing  is  done 
numerically  post  detection  after  A/D  conversion  of  the  detected  signal. 


Signal  Photons  Per  "Bit" 

Figure  8.  Basic  heterodyne  measurements  with  theory  and  simulation  results. 


Figure  8  plots  the  signal  to  noise  of  the  measurements  made  with  this  system  along  with  theoretical 
calculations  and  simulation  results.  All  signal  measurements  are  for  detected  photons  (not  incident 
photons).The  simulation  and  theory  are  in  complete  agreement.  The  actual  signal  to  noise 
measurements  suffer  from  added  thermal  and  amplifier  noise  and  are  about  17%  less  than  the 
theoretical  and  simulation  results.  When  the  input  lasers  are  turned  off  the  thermal  and  amplifier  noise 
is  easily  measured.  When  this  noise  is  subtracted  (in  quadrature)  from  the  measured  signal  only  the 
photon  shot  noise  should  remain.  The  resulting  S/N  is  only  about  5%  less  than  the  theoretical  and  within 
measurement  error  limits. 


2.3  Summary 

Simulations  of  90  and  180  degree  hybrids  were  carried  out  and  agree  well  with  theoretical  models. 
Measurement  of  180  degree  hybrid  signal  to  noise  ratios  agree  with  theoretical  and  simulations  within 
measurement  error.  In  the  coming  months  measurement  of  phase  noise  will  be  carried  out  on  the  180 
degree  hybrid  system  to  verify  simulations. 
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3  Long-Wave  Infrared  (LWIR)  Communications  Development 

(Mr.  David  Legare  and  Mr.  David  Hummel) 


3.1  Introduction 

The  focus  of  this  effort  is  on  the  evaluation  of  IR  spectral  region  in  the  8  to  14  micron  wavelength  range 
for  improved  free-space  communications  performance  over  legacy  short  wave  IR  (1.55  microns)  in 
adverse  atmospheric  conditions  including  aerosol  scattering  and  turbulence.  Improved  immunity  to 
aerosol  scattering  and  turbulence-induced  scintillation  at  longer  wavelengths  has  been  theoretically 
modelled  and  experimentally  verified  in  the  past.  However,  the  technology  required  to  engineer  a  long 
range,  high  data  rate,  low  SWAP  communications  system  in  the  LWIR  spectrum  has  only  very  recently 
become  an  emerging  possibility.  The  initial  phase  of  the  effort  described  below  entailed  a  set  of 
experiments  to  measure  relative  performance  over  a  short-range  link  through  a  homogenous 
atmosphere.  The  follow-on  phase,  currently  on-going,  will  conduct  long  duration  free-space  atmospheric 
propagation  measurement  and  performance  evaluation  under  operational  conditions  over  a  30  km  link 
between  AFRL-Rome  and  our  Stockbridge  test  site.  The  ultimate  program  goal  is  to  develop  and 
evaluate  a  prototype  LWIR  system  which  represents  a  solution  to  an  anticipated  real  operational 
requirement. 

3.2  Background 

Communications  in  the  optical  spectrum  offers  many  advantages  including  low  SWAP,  extremely  high 
bandwidth,  and  inherent  covertness  and  jam-resistance  due  to  spatial  confinement  of  the  radiated 
power  to  a  very  small  free-space  path  volume.  The  latter  also  helps  to  minimize  interference  and 
spectral  crowding  in  congested  environments.  Component  technology  covering  most  of  the  visible  and 
IR  spectrum  is  now  already  available,  or  at  least  reaching  the  point  of  foreseeable  maturity,  to  provide 
practical  implementations  of  compact,  very  high  speed  free-space  communications  links  to  support  a 
number  of  applications.  However,  the  physics  of  propagation  through  the  lower  atmosphere  will  always 
limit  the  utility  of  optical  communications  in  this  realm.  Therefore,  it  only  makes  sense  to  fully 
investigate  the  use  of  that  portion  of  the  spectrum  which  is  least  susceptible  to  these  detrimental 
effects. 

Fortunately,  as  can  be  seen  in  an  atmospheric  transmission  plot,  figure  1,  a  number  of  windows  t  across 
the  optical  and  IR  spectrum  exist  which  allow  for  low  loss  propagation  in  clear  air  conditions.  Since  the 
basic  molecular  composition  of  the  atmosphere  (minus  water  vapor)  will  always  remain  constant,  these 
are  the  only  usable  portions  of  this  spectrum  that  are  of  interest  for  all  but  very  short  range  links.  The 
essential  issue  then  becomes  energy  loss  via  scattering  from  aerosols,  mainly  comprised  of  water,  which 
are  distributed  along  the  propagation  path.  A  less  dominant,  but  still  very  significant  loss  mechanism  is 
turbulence-induced  scintillation.  This  will  be  investigated  more  thoroughly  in  the  follow-on  phase  of  the 
effort. 


Atmospheric  aerosols  may  be  defined  as  any  condensed  particle  of  water,  in  a  liquid  or  frozen  state 
either  falling  through  or  suspended  in  the  atmosphere.  As  such,  these  can  include  haze  and  fog,  with 
particle  sizes  ranging  from  a  fraction  of  a  micron  to  10  or  more  microns,  up  to  fine  mist,  rain,  snow,  sleet 
and  hail.  Scattering  is  greatest  when  the  wavelength  of  the  signal  is  on  the  same  order  as  the  diameter 
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of  the  aerosol  particle.  The  scattering  mechanism  is  effectively  the  absorption  and  radiative 
redistribution  of  the  signal  energy  in  a  pattern  which  is  peaked  in  the  direction  of  the  signal  entering  the 
aerosol.  When  the  particle  sizes  are  very  large  compared  to  the  wavelength,  this  scattering  effect  is 
minimal,  and  the  main  loss  mechanism  is  simple  absorption.  Thus,  in  rain  or  snow,  the  losses  at  both  the 
SWIR  and  LWIR  are  about  the  same.  This  is  apparent  in  Figure  10  showing  a  plot  of  signal  levels  taken 
during  a  snow  storm.  Here  it  can  be  seen  that  the  signal  levels  for  both  the  SWIR  and  LWIR  vary  about 
equally  with  a  change  in  precipitation  rate.  The  rate  of  snowfall  went  from  light  to  moderate,  then  back 
to  light  during  the  test  period. 


EM  SPECTRUM  CLEAR  AIR  TRANSMISSION  PLOT 


Figure  9.  Clear  air  transmission  plot. 
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Figure  10.  Attenuation  with  large  aerosols  (snow). 
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With  the  exception  of  clouds,  which  contain  a  very  dense  distribution  of  the  full  range  of  aerosol  sizes, 
haze  and  fog  are  predominantly  comprised  of  smaller  particles,  typically  less  than  a  few  microns  in 
diameter.  Since  it  is  extremely  common  for  at  least  some  fog  or  haze  to  be  present  along  a  portion  of  a 
signal  path,  it  is  reasonable  to  assume  that  operation  in  the  LWIR  spectrum  could  provide  for  improved 
information  throughput  and  link  availability  compared  to  legacy  SWIR  systems.  Even  on  apparently  clear 
days,  the  atmosphere  almost  always  contains  a  concentration  of  small  aerosols  that  can  significantly 
attenuate  an  optical  signal,  especially  over  long  distances.  So  it  can  be  seen  that  even  a  few  tenths  of  a 
dB/km  of  signal  loss  due  to  these  thinly  distributed  aerosols  can  have  an  appreciable  effect  over  long 
range  links  of  several  10's  of  kilometers.  Therefore,  although  it  makes  sense  to  employ  the  SWIR 
spectrum  for  space  and  high  altitude  applications,  there  are  some  real  advantages  to  utilizing  the  LWIR 
spectrum  when  any  portion  of  the  link  transgresses  the  lower  atmosphere. 

Operation  in  the  LWIR  spectrum  appears  to  provide  the  greatest  advantage  of  overcoming  adverse 
atmospheric  effects,  and  needs  to  be  thoroughly  investigated  for  employment  over  long-range 
communications  links.  Although  it  is  clear  that  no  optical  wavelength  will  operate  under  all  conditions, 
the  underlying  assumption  here  is  that  the  LWIR  will  perform  markedly  better  (on  the  order  of  at  least  a 
few  dB  per  km)  under  conditions  of  moderate  atmospheric  obscuration  (  haze,  short  patches  of  fog,  etc.) 
which  are  present  over  at  least  some  portion  of  a  long  path  length  a  high  percentage  of  the  time.  This 
makes  the  use  of  the  optical  spectrum  more  practical  by  allowing  for  a  much  higher  system  availability. 
The  exact  availability  numbers  are  not  yet  known,  as  this  analysis  is  part  of  the  reason  for  conducting 
this  program.  However,  it  is  assumed  that  they  would  show  a  significant  improvement  over  legacy 
optical  systems.  The  original  link  scenario  considered  for  this  program  assumed  current  ISR  downlink 
distances  of  around  10  to  30  Km.  Distances  of  100  miles  or  more  might  benefit  even  more  due  to  the 
fact  that,  due  to  earth  geometry,  most  of  the  path  would  occur  in  the  upper  atmosphere  where  aerosol 
concentration  is  low.  However,  even  a  small  difference  in  attenuation  over  a  very  long  distance  could 
make  or  break  the  link.  The  Mie  Scatter  models  and  the  literature  suggests  that  there  is  an  overall 
advantage  to  operating  in  the  LWIR  spectrum.  What  we  are  trying  to  do  here  is  prove  it  out  for  ourselves 
experimentally. 


The  major  objectives  of  the  program  are  to  quantify  the  relationship  between  atmospheric  parameters 
and  signal  attenuation,  and  to  make  a  comparison  between  the  performance  of  SWIR  wavelengths  (1.5 
microns)  and  LWIR  wavelengths  (10.6  microns)  under  similar  atmospheric  conditions.  Two  key 
atmospheric  attenuation  mechanisms  account  for  most  of  the  transmission  path  loss  in  the  IR  spectrum. 
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These  are  Mie  scattering  and  turbulence-induced  scintillation.  Mie  scattering  results  from  an  interaction 
between  the  IR  signal  and  aerosol  particles  such  as  those  microscopic  water  droplets  found  in  haze  and 
fog.  When  the  particle  diameter  is  less  than  or  approximately  equal  to  the  wavelength  of  the  impinging 
photons,  the  energy  is  essentially  redistributed  in  a  pattern  around  the  particle,  with  the  bulk  of  it 
travelling  forward  in  a  cone  of  small  angles  off  the  main  beam.  This  radiation  pattern  is  defined  as  a 
phase  function.  An  example  phase  function  plot  (for  some  given  particle)  is  shown  in  Figure  11.  Although 
most  of  the  signal  proceeds  along  its  original  path,  a  small  portion  is  sent  in  other  directions,  with  a 
resultant  loss  of  power  at  the  receiver. 


The  relationship  between  Mie  scattering  attenuation  and  aerosol  parameters  such  as  size  distribution 
and  concentration  is  fairly  complex.  The  extinction  coefficient  a,  given  in  equation  (3.1),  is  not  only  a 
function  of  wavelength,  but  of  physical  properties  of  the  individual  aerosol  particles.  It  is  obvious  from 
this  expression  that  direct  calculation  of  a  for  a  given  atmosphere  would  be  a  computationally  intensive 
task,  mainly  due  to  the  necessity  of  first  having  an  accurate  measure  of  the  aerosol  particle  size 

dN(r ) 

distribution  per  unit  volume  .  Thus,  the  real  problem  in  calculating  attenuation  using  this  formula 
is  in  obtaining  the  values  of  these  atmospheric  variables. 


rdr 


where  a  is  the  extinction  coefficient  in  Km'1,  A  is  the  wavelength  in  microns. 


f  2nr  ,3 

fdNl 

J>l 

A  ’n 

v  A  y 

dr  J 

(3.1) 


dN(r ) 


is  the  particle  size 


distribution  per  unit  volume  (in  cm'4),  n'  is  the  real  part  of  the  refractive  index  of  the  aerosol,  r,  is  the 
radius  of  a  given  particle  size  in  cm,  and  Qd  ( —  ,n'J  is  the  scattering  cross-section  for  a  given  aerosol 
type. 


It  is  extremely  difficult  to  measure  aerosol  particle  size  distributions  and  concentrations  at  any  given 
point,  and  especially  over  a  long  distance  where  their  values  can  vary  dramatically  in  both  time  and 
space.  The  goal  of  the  current  program  is  not  an  attempt  to  make  such  measurements  or  to  develop  a 
scientific  understanding  of  the  physics  of  these  scattering  mechanisms.  Rather,  we  correlate  link 
performance  with  an  easily  obtainable  atmospheric  parameter  that  directly  relates  to  the  bulk  effect  of 
these  aerosols  along  the  transmission  path.  The  Kruse  model,  discussed  in  the  next  section,  provides  an 
experimentally-verified  expression  for  Mie  scattering  based  on  meteorological  visibility  V.  As  such,  the 
extinction  coefficient  a  can  be  expressed  for  any  given  wavelength  A.  as  a  function  of  a  single  variable  V. 


Although  V  is  often  determined  crudely  as  the  distance  at  which  a  human  observer  can  discern  a  certain 
image  (typically  through  haze  or  fog),  it  is  quantitatively  defined  as  the  distance  travelled  by  an  optical 
beam  before  it  is  reduced  to  5%  of  its  original  value.  The  current  program  will  correlate  experimentally- 
derived  attenuation  of  simultaneously  transmitted  SWIR  and  LWIR  signals  over  a  folded  path  at  near 
ground  level  with  observed  measures  of  atmospheric  visibility.  Visibility  estimates  used  to  date  have 
been  based  on  a  combination  of  on-line  reported  values  taken  from  the  local  area,  as  well  as  direct 
human  observation.  These  results  will  also  be  correlated  with  values  derived  from  the  Kruse  model.  This 
will  provide  an  additional  level  of  confidence  in  the  validity  of  the  experiment,  assuming  the  empirical 
and  theoretical  results  are  in  agreement.  It  will  also  serve  to  verify  the  utility  of  the  Kruse  model  for 
extension  into  the  IR  portion  of  the  spectrum,  since  it  was  primarily  derived  from  data  collected  in  the 
visible  portion  of  the  optical  spectrum.  Note  that  although  the  results  reported  here  were  made  based 
on  the  methods  mentioned  previously,  we  have  since  obtained  an  actual  visibility  meter.  Figure  12, 
which  will  be  available  for  use  in  future  data  collection.  This  device  works  on  the  principle  of  off-axis 
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scattering  (particle  phase  function)  between  an  optical  emitter  and  detector  separated  by  a  few  feet. 
The  values  provided  would  be  accurate  for  our  purposes  assuming  that  the  visibility  measured  over  this 
short  distance  is  uniform  over  the  entire  transmission  path.  It  is  safe  to  assume  that  this  is  the  case  for 
our  .78Km  link. 

Figures  13-16  are  images  of  visibility  used  to  estimate  V  for  data  taken  during  the  experiment.  Figures  13 
and  14  are  a  view  looking  West  from  the  shelter  housing  the  experiment  transceiver.  The  distance  to  the 
far  horizon  in  this  view  is  about  1  Km.  Figure  13  is  the  view  under  clear  sky  conditions.  In  Figure  14,  it 
can  be  seen  that  the  most  distant  tree  line  is  mostly  obscured.  This  image  therefore  corresponds  to  a 
visibility  value  V  of  about  1  Km.  A  similar  comparison  can  be  seen  in  the  two  views  shown  in  figures  15 
and  16.  These  are  views  looking  East  from  the  same  shelter.  With  figure  15  as  the  clear  sky  view,  and 
figure  16  likewise  reveals  an  obscured  tree  line  at  the  horizon,  a  distance  of  about  2  Km  from  the 
shelter.  Other  values  of  V  for  the  purposes  of  the  experiment  were  based  on  visual  estimates  using 
these  two  as  a  reference. 


It 
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Figure  14.  Overcast  sky  visibility 
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Figure  15.  Clear  sky  visibility 
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Figure  16.  Overcast  sky  visibility 
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3.3  Tasking 


This  report  describes  the  initial  (baseline)  effort  to  conduct  atmospheric  propagation  testing  to  evaluate 
the  comparative  performance  of  the  LWIR  spectrum  (at  around  10.6  microns)  and  the  SWIR  spectrum 
(at  around  1.5  microns)  under  a  range  of  atmospheric  conditions  including  rain,  snow,  fog,  and  haze.  The 
test-bed  was  established  over  a  horizontal  short-range  folded  path  comprised  of  a  retro-reflector  placed 
1300  feet  away  from  the  transceiver  hardware  to  provide  for  a  total  path  length  of  about  .78  Km.  The 
homogeneity  of  the  atmosphere  over  this  short  range  provided  for  a  more  accurate  correlation  between 
the  observed  visibility  parameter  V  and  the  measured  signal  values  propagating  through  this  well- 
characterized  atmospheric  volume.  A  diagram  of  the  transceiver  hardware  is  shown  below  in  Figure  17. 


TransmitAperture 


Figure  17.  Diagram  of  optical  transceiver  assembly 


3.4  Theoretical  Model 

A  link  budget  model  was  developed  to  theoretically  characterize  the  propagation  for  both  the  LWIR  and 
SWIR  signals.  Equation  (3.2)  is  the  power  received  from  the  transmitter.  The  model  can  essentially  be 
broken  down  into  three  major  loss  factors:  1)  Free  space  path  loss,  2)  Losses  due  to  aerosol  scattering, 
and  3)  System  losses.  Variables  entering  into  the  link  budget  calculation  are  listed  in  Tables  2  and  3  The 
scatter  loss  portion  was  based  on  the  Kruse  model,  which  calculates  the  extinction  coefficient,  a,  as 
shown  in  equation  (3.3). 

Figure  18  is  a  plot  of  the  of  the  link  budget  calculated  for  the  .78  Km  test  link  used  in  the  experiment. 
Note  that  the  system  parameters  entered  for  the  model  were  tailored  so  that  both  the  SWIR  and  LWIR 
signals  would  have  the  same  received  power  at  maximum  visibility,  where  losses  due  to  scattering  were 
nil.  Thus,  since  equal  aperture  sizes  and  transmitted  power  levels  at  both  wavelengths  would  minimize 
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geometric  losses  (free  space  path  loss)  for  the  SWIR  signal,  these  values  were  reduced  for  the  SWIR  so 
that  only  the  comparative  losses  due  to  scattering  would  be  apparent. 
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Table  2  Extinction  Exponents  vs  Visibility 


Visibility  V 

V  >  50 

6  <  V  <  50 

1  <  V  <  6 

0.5  <  V  <  1 

V  <  0.5 

q 

1.6 

1.3 

0.16  V+ 0.34 

< 

i 

Ln 

0 

Table  3  Link  Budget  Parameters 


D  =  Rx  Aperture  Diameter  (m) 

dt  =  Tx  Aperture  Diameter  (m) 

L=  Path  Length  (m) 

X  =  Wavelength  (m) 

T,=  Tx  System  Loss 

Tr  =  Rx  System  Loss 

Visibility 

A 

Visibility 

A 

(km) 

1550nm 

(km) 

10600 nm 

0.8 

-11.8851 

0.8 

-5.75859 

0.9 

-10.5765 

0.9 

-5.11432 

1 

-9.53184 

1 

-4.60615 

1.1 

-8.67914 

1.1 

-4.19664 

1.2 

-7.97037 

1.2 

-3.86088 

1.3 

-7.37229 

1.3 

-3.5816 

1.4 

-6.86116 

1.4 

-3.3465 

1.5 

-6.41956 

1.5 
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Figure  18.  Comparative  link  budgets:  SWIR  vs  LWIR 
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3.5  Method  of  Data  Collection 

The  purpose  of  this  first  phase  of  the  LWIR  in-house  program  was  to  accurately  characterize 
atmospheric  propagation  loss  due  to  aerosol  scattering  simultaneously  at  both  the  SWIR  (1.55  p)  and 
the  LW  IR  (10.6  p).  The  desired  result  was  to  quantify  both  the  relative  performance  of  the  two 
wavelengths,  and  the  absolute  performance  of  the  longer  IR  wavelengths  for  the  purpose  of  ultimately 
designing  an  "all  weather"  optical  communications  system  that  could  be  deployed  in  an  operational 
environment. 


The  aerosol  scatter  loss  measurements  will  were  made  over  a  folded  path  using  a  5"  diameter  first- 
surface,  gold-coated  retro-reflector  about  1300  feet  from  the  transmitter/receiver  system  for  a  total 
distance  of  about  .78  Km.  The  transmitter/receiver  assembly  consisted  of  a  3'x3'  optical  breadboard 
containing  the  optical  and  electro-optic  hardware,  along  with  the  control  and  measurement  electronics 
and  software  required  to  generate  both  of  the  test  wavelengths.  The  test  system  was  mounted  on  a 
precision  pointing  gimbal  and  co-aligned  with  a  high  resolution  telephoto  camera  so  that  it  could  be 
precisely  bore-sighted  on  the  retro-reflector  at  the  far  end  of  the  propagation  path.  Measurements 
were  taken  at  all  possible  opportunities  under  a  wide  range  of  atmospheric  conditions  to  include  clear 
air  (maximum  visibility)  and  fog/haze,  rain  and  snow.  The  test  range  was  located  at  the  AFRL/Rome  site 
and  conducted  along  a  horizontal  ground  level  path.  The  transmitter/ receiver  system  was  operated 
from  inside  an  environmentally-controlled  military  field  shelter  (through  an  open  door),  so  that  it  could 
be  conveniently  protected  while  experiments  were  conducted  under  adverse  weather  conditions.  The 
shelter  was  located  in  the  Satcom  test  area  on  the  South  side  of  Building  3  at  AFRL/RI,  in  Rome,  NY.  The 
other  end  of  the  folded  path  was  located  at  the  far  southwest  end  of  the  building  3  parking  lot,  about 
1300  feet  away.  The  retro-reflector  was  positioned  on  a  20'  tower  attached  to  our  FIMMWV  at  this 
point. 


Both  the  SWIR  (1.55  p)  and  LWIR  (10.6  p)  signals  were  generated  simultaneously  at  a  5  KHz  pulse  rate 
and  emitted  along  a  parallel  path  from  the  transmit  side  of  the  test  system.  The  signals  returned  from 
the  retro-reflector  were  collected  into  a  5"  aperture  telescope  which  was  mounted  next  to,  and  aligned 
with  the  transmit  aperture.  The  output  of  the  receive  telescope  proceeds  through  a  beam  splitter  to 
separate  the  two  wavelengths  into  separate  paths  for  subsequent  beam  shaping,  optical  filtering  and 
detection.  The  detector  outputs  were  fed  to  separate  channels  of  an  oscilloscope  which  was 
programmed  to  average  the  peak  amplitudes  of  each  pulse  (in  millivolts)  in  the  train  with  about  30 
seconds  of  averaged  pulses  comprising  each  data  point.  A  number  of  these  data  points,  taken  over  a 
typical  period  of  several  minutes  to  an  hour  during  which  the  atmospheric  conditions  (ie  fog/haze 
density)  remained  relatively  constant,  were  then  averaged  to  provide  a  single  value  for  received  power 
corresponding  to  the  observed  value  of  constant  atmospheric  visibility  during  that  time  period.  These 
values  represent  the  expected  received  power  levels  for  that  given  atmospheric  condition,  and 
subsequently  were  used  to  calculate  an  aerosol  extinction  loss  for  each  wavelength  in  dB  as: 


LExt  =10Log 


PR 

PT 


(  PR 

-10  Log1  cal 


V  PTcal  J 


0-4) 


Where:  PT,  PR  =  Transmit,  Receive  Power  for  a  given  test  during  non-clear  air  conditions,  and  PTcal, 
PRcal  =  Transmit,  Receive  power  measured  during  a  clear  air  (maximum  visibility;  no  aerosol  scattering 
loss)  test.  The  latter  was  used  as  a  calibration  value  which  essentially  provides  a  measure  of  the 
combined  transmitter/receiver  system  loss  and  free-space  path  loss.  Therefore,  the  second  term  of  LExt 
represents  the  combined  fixed  test  system  loss  plus  the  total  clear  air  path  loss,  both  of  which  should 
never  change.  Thus,  this  serves  as  a  reference  from  which  to  derive  the  actual  value  of  just  the  scatter 
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extinction  loss  for  any  given  atmospheric  condition.  Since  the  atmospheric  path  length  is  fixed  at  .78  Km, 
the  extinction  loss  for  the  given  atmospheric  test  can  be  easily  expressed  in  standard  form  as: 

Aerosol  Extinction  Loss  =  LExt/.78  dB/Km  (3.5) 

Note  that  PT  and  PTcal  values  used  in  the  calculation  were  measured  as  small  signals  scattered  off  a 
beam  splitter  in  the  transmitter  chain  in  real  time,  and  therefore  were  only  a  small  unknown  fixed 
fraction  of  the  total  transmitted  power.  This  still  makes  for  a  valid  calculation  since  the  ratio  of  this 
measured  power  to  actual  total  transmitted  power  is  unimportant,  as  long  as  the  ratio  remains  constant 
throughout  the  test.  Thus,  the  laser  sources  were  controlled  so  that  these  PT  and  PTcal  values  were 
kept  constant  throughout  all  of  the  experiments. 

Also  note  that  the  signal  transmit  power  levels  for  both  lasers  were  adjusted  so  that  the  receive  power 
for  each  wavelength  was  close  to  equal  under  clear  air  conditions.  Although  this  wasn't  necessary  to 
make  the  actual  loss  measurements  (since  the  loss  at  each  wavelength  is  calculated  as  a  ratio  between 
its  own  transmit  and  receive  value),  it  allows  for  a  more  direct  comparison  of  the  data  for  each 
wavelength. 


3.6  Data  Collected 

Due  to  initial  problems  with  the  experiment  hardware  and  the  inability  to  continue  the  data  collection 
effort  beyond  August  3013,  it  was  only  possible  to  obtain  a  few  good  data  sets  during  this  first  phase  of 
the  effort.  The  first  of  these  shown  in  Table  4  below  were  taken  over  a  series  of  three  clear  air  days 
where  the  atmospheric  visibility  was  known  to  be  at  a  maximum.  The  results  from  each  of  these  days 
were  averaged  together  to  provide  a  best  value  for  PRCal  at  each  wavelength.  The  reference  values, 
shown  as  287  mv  (for  X  =  1550  nm)  and  307  mv  (for  X  =  10600nm)  are  used  hereinafter  to  calculate  the 
extinction  loss  LExt  for  data  taken  at  other  (lower)  values  of  atmospheric  visibility  shown  in  Tables  5, 
3.15,  and  3.6.  It  can  be  readily  seen  that  they  are  well  in  agreement  with  the  results  predicted  by  the 
Kruse  model.  What  is  particularly  interesting  is  the  very  close  correlation  between  the  comparative  path 
loss  values  (A  dB/Km)  derived  from  the  empirical  data  and  what  is  predicted  by  the  model  in  figure  15. 
Thus,  though  less  than  the  originally  desired  amount  of  data  was  taken,  it  agrees  extremely  well  with  the 
predicted  results.  This  leads  to  a  very  high  level  of  confidence  in  the  Kruse  model  for  use  as  a  predictor 
of  atmospheric  aerosol  extinction  for  the  entire  IR  spectral  band  of  interest,  at  least  for  visibility  ranges 
of  around  one  Km  to  clear  air  conditions.  This  is  particularly  encouraging  since  this  simple  model  is 
based  solely  on  the  parameter  V  (atmospheric  visibility),  as  long  as  the  aerosols  are  essentially 
microscopic  in  size  (not  rain  or  snow).  Note  from  the  previous  background  discussion  (and  figure  10) 
that  the  primary  effect  of  these  larger  particles  is  absorption  (not  scattering),  and  that  they  have  an 
essentially  equal  effect  on  all  wavelengths  of  interest. 

Although  we  are  now  in  the  second  phase  of  the  effort  which  involves  data  collection  over  the  30  Km 
Griffiss-Stockbridge  link,  we  plan  to  re-establish  and  collect  additional  data  on  the  .78  Km  link  to  further 
characterize  and  validate  these  assumptions,  and  perhaps  provide  enhancements  to  the  Kruse  model 
based  on  observed  propagation  results  in  the  IR  spectrum. 
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Table  4  Clear  Air  Data  (for  measurement  calibration) 


.1  T 


25  April 


256.08 

257.73 

237.82 

241.34 

260.37 

253.11 
261.14 
249.04 
257.84 
239.58 
251.02 
260.48 

230.12 
248.16 
245.19 

239.47 
242.55 

259.38 

249.47 


285.1 

314.8 

326.3 

322.3 

328.3 

297.1 

292.3 

308.9 

276.6 

302.7 
295.5 

304.7 

305.9 

309.2 


248.27 

250.58 

267.96 

289.96 
290.29 
284.35 
323.73 


301.18 

266.75 
267.85 
267.52 
330.44 

316.58 
325.27 

357.83 

334.84 

363.55 
339.68 
296.01 

310.75 
299.42 
309.54 

289.85 
330.66 
318.01 

310.86 

305.58 

307.56 
312.84 
283.36 


10600  nmRx 

250.1 
231.3 

234.9 

244.2 

240.5 

234.2 

246.9 

249.2 

255.3 
275 

282.2 

297.3 

286.3 

296.9 

295.5 

316.9 

337.7 

329.6 

323.1 

342.9 

331.7 

310.2 

317.8 

357.4 

360.4 

293.4 

335.6 

344.4 

342.8 

324.8 

312.1 

332.1 

352.9 

374.1 

304.9 

364.2 

323.8 

296.4 
328.1 

322.5 

290.7 


1550  nm  R> 

358.3 

333.8 

362.4 

332.4 

342.9 
285.1 
351.3 


279.3 

299.4 

275.5 

263.5 

289.8 

286.4 

274.2 

268.4 

270.6 

262.6 

324.9 

330.6 

371.7 

398.7 

377.3 

294.7 
335.6 

339.9 


10600  nmF 

327.8 

320.9 

312.9 
311.1 
270 
247 


319.8 

283.6 

306.1 

307.6 

325.1 

275.1 

290.3 

294.3 
319.8 

303.6 


Reference  Clear  Air 


Table  5  Data  at  3  Km  Visibility  (5/10/13) 


Table  6  Data  at  1  Km  Visibility  (8/30/13) 


10  May  2013 


Time 

1.55  Rk 

10.6  Rk 

16 :1b 

16/97 

262 

16:19 

185.9 

286.7 

16:21 

197.34 

237.S 

16:22 

189.53 

24  7.  S 

1fi:?4 

173.36 

755  fi 

16:26 

157.19 

270  7 

16:28 

167.2 

285.2 

16:29 

187-27 

789  3 

Average 

178 

209 

Cuinpur  utivv  Path  Loss  A((J8/Kni) 
Niwsurod:  1.92 

Krusn  Mndftl:  1.9? 


Ceneral  Atmospheric  Conditions 

Liyjhl  ruin 
Ligjil  hfi» 

Est  Visibility:  3  Km 

libOinn 

Path  LossUuo  to  Scattering 
Measured:  {d  F^/kmrJ  -?.fp[i 
Kruse  Model  [Ulsr  3  Km)  ?.4? 

IflfiOflnm 

Path  I  oss  One  to  Scattering 
Measured:  [dB/Km]  =  -.74 
Kruse  MnrieE  [Vi*:  3  Km)  =  -.SO 


30  Aug  2013 


Time 

1.55  Rk 

10.6  Rk 

8:55 

19.811 

106 

8:57 

18.414 

104.6 

9:00 

15.246 

104.2 

9:03 

31.559 

111.8 

9:04 

34.903 

116 

9:06 

45.353 

122.4 

9:07 

49.852 

126.1 

9:08 

44.913 

126.9 

9:09 

40.546 

119.6 

9:10 

39.721 

117.7 

34.0318 

115.53 

Average 

34 

116 

General  Atmospheric  Conditions 

light  milt 
Moderate  Fog 
1st  Visibility:  1.0  Km 

ISSOnm 

Path  l  o»s  Oue  to  Scattering 
Measured:  (dB/Vm)  *  -11.87 
Kruse  Model  (Viv  1.0  Km)  • -10.13 

lOGOQnm 

Path  Loss  Due  to  Scattering 
Measured:  (d8/Km)  » -5.42 
Kruse  Model  (Vis:  1.0  Km)  -  3.9 


Comparative  Path  loss  A  (dB/Km) 
Measured:  6.45 
Kruse  Model:  6.23 


Table  7  Data  at  3  Km  Visibility  (5/23/13) 


23  May  2013 

GepbLT  al  AtmeiipliEriC  Cortdlli£w*£. 

Time 

1.55  Rk 

10.6  fts 

Light  rain 

9  01 

13972 

?5M 

Ls.|fni 

4  If;.  -Z  fern 

907 

107.09 

255-9 

9:04 

177.37 

200-5 

ISSOnm 

Pjlfth  1  ARC  IW  ' 

9:0ft 

156.31 

t-.-i  Crjiftsrina 

9:30 

1*9.64 

245? 

J.5i '.l'.ijrt'tl :  IdR/kml  i  .2.72 

9:22 

164.78 

233-2 

Kruw  Model  I'Vk:  .1  Kmi 

a  -2-4? 

9:24 

151,76 

250-5 

93? 

149.71 

274-3 

lOOOOnm 

9:38 

131.14 

275-0 

Path  Loss  Pue  to  Scatteri  ng 

9:30 

15961 

260-4 

MrdSwed:  (dB/Kinl  ..57 

9:32 

173.91 

2654 

9:37 

171.38 

2B4.3 

9:38 

166.1 

2E5.3 

9:40 

170.5 

2E2.7 

9:46 

187.55 

264. 2 

10:01 

194.48 

293 

10:10 

170.5 

265.7 

10:35 

236.61 

326.2 

10:17 

325  72 

348.5 

10:54 

191  02 

3344 

Average 

170 

277 

Lcmp&rative  Path  LOS*  A  (dB/tort) 

Mi.MLHeR'd:  2.1$ 

Krase  Model:  1.92 
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3.7  Analysis/Conclusions 

3.7.1  SWIR 


There  are  many  considerations  in  designing  an  atmospheric  optical  communications  link.  An  obvious 
advantage  of  legacy  SWIR  (1.5  p)  communications  is  about  a  16dB  margin  in  geometric  space  path  loss 
(due  to  higher  gain  resulting  from  a  higher  optical  frequency)  over  the  LWIR  (at  10.6  p)  for  the  same  size 
transmit/ receive  aperture.  However,  the  resulting  narrower  beam-width  also  makes  pointing  more 
difficult  due  the  required  increase  in  pointing  accuracy,  and  the  random  beam  wander  resulting  from 
turbulence-induced  scintillation.  These  scintillation  effects  also  introduce  jitter  on  the  received  signal  in 
the  detector  focal  plane,  which  can  further  result  in  an  average  signal  degradation,  and  lowered  signal 
to  noise  ratio.  Longer  optical  wavelengths  are  much  less  susceptible  to  all  of  these  scintillation  effects, 
and  thus  have  a  countering  advantage  which  may  overcome  at  least  some,  if  not  all  of  the  obvious 
geometric  path  loss  advantages  of  an  SWIR  system.  These  effects  will  be  investigated  further  in  the  next 
phase  of  the  effort  using  data  collected  over  the  30  Km  Griffiss-Stockbridge  link.  The  purpose  of  the 
present  analysis  is  to  quantify  the  effects  of  aerosol  scattering.  The  data  collected  showed  a  good 
correlation  with  the  Kruse  model  for  visibilities  ranging  from  around  1  Km  to  clear  air  conditions.  This 
model  will  now  be  used  to  predict  the  comparative  performance  of  the  SWIR  and  LWIR  wavelengths 
over  a  hypothetical  50  Km  link. 

Figure  19  shows  the  effects  of  scattering  alone.  In  this  case,  the  transmit  and  receive  aperture  sizes  were 
selected  such  that  the  link  margins  would  be  equal  at  very  high  visibility,  where  the  scattering  losses 
were  nil.  For  comparison,  in  figure  20,  the  model  was  run  with  equal  aperture  sizes  for  both 
wavelengths  to  more  fully  illustrate  all  components  of  the  link  budget.  Here,  the  16dB  SWIR  margin 
becomes  apparent  as  the  plot  approaches  100%  percent  atmospheric  visibility,  and  the  scattering  loss 
becomes  negligible.  Figure  21  further  illustrates  a  case  nearly  optimized  for  the  LWIR  in  which  the  10.6  p 
wavelength  operates  with  an  aperture  of  twice  the  diameter  of  the  SWIR  to  close  the  gap  in  the 
geometric  path  loss  advantage.  Though  doubling  may  at  first  seem  like  a  considerable  change,  it  is  seen 
that  it  is  only  an  increase  from  5  cm  to  10  cm  (about  2"  to  4"). 

All  else  being  equal,  it  is  clear  from  figure  21.  that  the  LWIR  has  a  performance  advantage  where  the 
average  visibility  over  the  link  is  less  than  about  15  Km.  We  know  from  our  own  experience  that  it  is  very 
difficult  to  differentiate  visibility  levels  beyond  about  6  to  7  kilometers  from  clear  air  conditions  purely 
by  sight.  Therefore,  it  would  be  equally  difficult  to  estimate  the  percentage  of  time  that  the  link  would 
experience  average  visibility  conditions  exceeding  this  15  Km  threshold  without  taking  actual  long  term 
data  measurements.  This  will  be  done  in  the  next  phase  of  the  effort. 

Thus,  a  major  goal  in  this  next  phase  will  be  to  quantify  the  atmospheric  statistics  that  would  relate  to 
the  availability  of  an  operation  link  over  a  significant  path  length  of  at  least  a  few  tens  of  kilometers.  A 
distinguishing  feature  of  longer  links  is  the  variability  of  the  atmospheric  parameters,  including  visibility, 
along  the  path.  Note  that  the  data  collected  in  the  present  phase  was  over  a  relatively  short  (.78  km) 
folded  path  over  which  the  atmospheric  conditions  could  be  assumed  to  be  homogenous.  Ultimately,  it 
would  be  desirable  to  express  the  varying  visibility  over  a  much  longer  path  as  a  single  value  of 
equivalent  (or  average)  visibility  which  could  be  correlated  with  the  Kruse  model.  Values  for  link 
availability  could  then  be  derived  as  a  function  of  the  atmospheric  statistics  along  the  path. 

Data  collection  will  take  place  over  the  wide  range  of  seasonal  atmospheric  conditions  that  normally 
occur  in  the  Central  New  York  area.  The  results  provided  by  the  Griffiss-Stockbridge  test-bed  could 
therefore  support  a  first-order  prediction  of  SWIR  and  LWIR  absolute  and  comparative  performance  for 
a  variety  of  anticipated  user  applications  in  a  wide  range  of  global  operating  environments. 
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50  Km  Link  Budget 

Atmospheric  Visibility  (km) 

2  4  6  8  10  11  11.7  12.2  13  15  15.5  18  22  30 


-200 


- SWIR  1550  nm 

- IW1R  10600  nm 


Visibility  (km)  LWIR  Prec/Ptx  (dB)  SWIR  Prec/Ptx  (dB) 

2  -91.42422  -248.57 

3  -55.46466  -153.378 

4  -41.85772  -108.335 

5  -35.87258  -83.039 

6  -33.01356  -67.3964 

7  -31.57529  -57.1101 

8  -30.82572  -50.0529 

9  -30.42495  -45.0592 

10  -30.20649  -41.4418 

10.5  -30.13728  -40.0068 

11  -30.08561  -38.7729 

11.2  -30.06879  -38.3267 

11.7  -30.03431  -37.3163 

12  -30.0179  -36.7748 

12.2  -30.00841  -36.438 

12.5  -29.99606  -35.9663 

13  -29.9796  -35.2607 

14  -29.95775  -34.1018 

15  -29.94521  -33.2074 

15.2  -29.94342  -33.0541 

15.5  -29.94109  -32.8381 

16  -29.93796  -32.512 

18  -29.93129  -31.5403 

20  -29.92901  -30.9339 

22  -29.92822  -  30.5494 
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Figure  19-  Scatter  Loss  Only 
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Figure  20  -  Equal  Aperture  Size 
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Figure  21  -  Larger  LWIR  Aperture 
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3.7.2  MWIR  to  LWIR  Band 


The  simple  Kruse  model  has  been  demonstrated  to  be  a  reasonably  accurate  predictor  of  scattering 
performance  in  the  LWIR  band  (around  10.6  microns)  and  the  SWIR  band  (around  1.5  microns).  It  would 
thus  be  safe  to  assume  that  this  model  could  be  used  to  estimate  comparative  performance  of 
wavelengths  between  these  extremes,  assuming  that  they  are  also  within  a  clear-air  transmission 
window.  Such  conditions  exist  for  the  spectral  regions  of  approximately  3  to  4.1,  4.3  to  5,  5  to  5.1,  and 
7.7  to  8.5  microns. 

Besides  the  obvious  advantage  of  higher  signal  gain  for  a  given  size  optical  aperture,  operation  at  these 
intermediate  wavelengths  also  has  a  potential  advantage  in  the  implementation  of  the  receiver.  To  date, 
the  most  promising  detector  technology  for  IR  detection  in  the  mid  to  far  IR  band  appears  to  be  the 
HgCdTe  eAPD  (Mercury  Cadmium  Telluride  electron  Avalanche  Photodiode).  This  technology  has  made 
significant  advances  in  the  past  several  years.  These  devices  in  general  seem  to  provide  the  best  overall 
speed  and  sensitivity  (D*)  over  all  known  commercially  available  detector  materials  and  configurations. 
However,  with  the  present  state-of-the-art,  the  advantage  is  significantly  less  at  wavelengths  beyond 
about  8  microns.  Thus,  at  least  in  the  near  term,  it  would  be  advantageous  to  consider  operation  at 
these  Mid  IR  wavelengths.  It  should  be  noted,  however,  that  we  also  have  developed  an  approved  SBIR 
topic  for  the  investigation  and  demonstration  of  advanced  LWIR  detector  technologies,  to  be  released  in 
FY15. 

An  analysis  of  some  of  these  intermediate  wavelengths  using  the  Kruse  model  also  indicates  that  the 
greater  part  of  the  propagation  gain  (reduced  scattering  through  atmospheric  aerosols)  demonstrated 
at  10  microns  can  be  achieved  with  wavelengths  in  the  Mid  IR  band,  as  shown  in  figure  22.  QCL  laser 
technology  is  also  readily  available  at  these  wavelengths,  making  the  construction  of  a  deployable 
system  a  reasonable  near-term  goal.  The  addition  of  long  range  propagation  testing  of  selected  Mid  IR 
wavelengths  over  the  Griffiss-Stockbridge  test  bed  is  tentatively  planned  for  FY15  or  FY16. 
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Figure  22.  Multiple  Band  Plot 
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4  Toward  an  Integrated  RF/FSO  Link 

(Mr.  John  Malowicki  and  Mr.  David  Hummel) 

4.1  Introduction 

A  communication  link  for  the  purposes  of  system  evaluation  and  experimental  testing  has  been 
established  by  AFRL/RI  between  its  Rome  Research  Site  (RRS)  and  Stockbridge  Test  Site  (STS).  It  includes 
an  optical  link  based  on  an  infrared  (1550  nm)  laser  carrier,  and  an  VW  RF  link  in  the  71-76  and  81-86 
GHz  bands.  In  the  near  future  a  Long  wave  IR  capability  will  be  added;  see  Section  3.  This  multi-band 
capability  will  serve  to  collect  long  term  link  performance  data  in  all  weather  conditions.  The  analysis  of 
this  data  will  provide  a  better  understanding  of  link  availability  and  throughput  which  will  help  to  set 
realistic  expectations  of  the  performance  of  future  fielded  systems. 

4.2  System  Overview 

The  RRS  node  is  mounted  on  an  80'  tower,  immediately  Southeast  of  Building  3  (N  43.22066°,  W 
75.40755°),  at  an  altitude  of  approximately  550  ft  above  sea  level.  The  Stockbridge  node  is  temporarily 
mounted  on  scaffolding  (N  43.02933°,  W  75.64585°).  By  the  end  of  2014  a  70'  tower  will  be  installed  at 
this  location  and  the  units  will  be  mounted  on  the  tower.  Operating  altitude  is  approximately  1350  - 
1380  ft  above  sea  level.  Figure  23  shows  the  two  locations  and  line  of  site  in  a  map  view.  Figures  24 
and  25  show  aerial  views  of  the  RRS  and  Stockbridge  sites,  respectively.  The  ground  elevation  profile  is 
shown  in  figure  26.  Accounting  for  tower  heights,  the  elevation  angle  from  the  RRS  node  to  the  STS 
node  (at  the  highest  tower  position)  is  0.56  degrees. 


Figure  23.  Map  of  Line  of  Sight  Optical  Path 
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Figure  24.  Rome  Site 


Figure  25.  Stockbridge  Site 
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Figure  26.  Elevation  Profile 


4.3  Parameters  of  the  Optical  System 

The  optical  terminal  is  manufactured  by  AOptix  and  is  similar  to  the  units  developed  under  the 
DARPA  FOENEX  program. 

•  Beacon  Laser  Diode  Array 
o  Wavelength:  850  nm 

o  Optical  power:  10  mW  (reduced  to  1  mW  with  ND  filter) 
o  Beam  Divergence:  0.2  deg  full  angle 
o  Path  distance:  28  km 
o  Beam  size  at  receiver:  ~100  m  diameter 

•  Communications  Laser 

o  Wavelength:  1545.31  and  1555.74  nm 

o  Optical  power:  Maximum  1.5  W  output  (but  restricted  to  1  W  for  increased  safety) 
o  Transmit  Optical  Aperture:  10  cm  diameter,  circular 
o  Beam  Divergence:  ~16  prad  (diffraction-limited) 
o  Path  distance:  28  km 
o  Beam  size  at  receiver:  ~45  cm  diameter 

4.4  Parameters  of  the  RF  System 

The  VW  system  is  commercially  available  from  BridgeWave  Communications  and  is  a  model  AR80 
with  24"  dish. 
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•  Frequency:  71  -  76  GHz  and  81  -  86  GHz 

•  Link  Budget: 

o  Gigabit  Ethernet  Mode:  186  dB  @  10 12  B.E.R. 

•  Antenna  Parameters 

o  Size:  24"  (60  cm)  Parabolic 
o  Main  Lobe  Beam  width:  0.4  degrees 
o  Peak  Gain:  23  dB 

4.5  Link  Budget  -  Optical  System 

A  first  approximation  link  budget  can  be  calculated  in  eq.  4.1  using  simple  geometrical  and  atmospheric 
attenuation  assumptions.  Figure  27  illustrates  the  geometry  of  the  transmitted  beam  and  the  collecting 
aperture  at  the  receiver.  Using  the  diameter  of  the  transmit  optics  and  assuming  a  divergence  D,  a  ratio 
of  the  area  of  the  collection  optics  to  the  beam  size  at  the  receiver  can  be  calculated.  This  ratio 
multiplied  by  the  power  transmitted  and  a  loss  factor  for  atmospheric  attenuation  represented  by  a, 
gives  an  estimate  of  the  received  optical  power.  The  value  of  the  atmospheric  attenuation  can  be 
estimated  from  Table  1  based  on  the  optical  visibility  of  the  atmospheric  path. 


received 


—  p  * 

1  transmitted 


d ; 


(4,  +  (D  *  R)): 


*K)(-“**/10) 


(4.1) 
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Table  8 


HfearAfr  Amount 

Condition  Precipitation  mm/hr  Visibility  dB  Loss/km 

Om,  50  m  -271 .65 


Dense  Fog. 

Thick  fog 
Moderate  log 
Liglu  Fog 

snow- 

snow 

Thin  Log 

snow- 

Haze 

snow 

Light  haze 

snow- 

Clear 

snow 

Very  Clear 

Cloudburst 

100 

Heavy  rai  n 

25 

Medium  rain 

12.5 

Light  rain 

25 

Drizzle 

0.25 

200  m 

-59.57 

500  m 

-20.99 

770  m 

-12,65 

1  km 

-9,26 

1 .9  km 

-4,22 

2  km 

-3,96 

2-8  km 

-2.58 

4  km 

-1,62 

5,9  km 

-0.96 

10  km 

-0.44 

18,1  km 

-0-24 

20  km 

-0.22 

25  km 

-11.19 

50  km 

-0,06 

Figure  28.  Optical  terminal  and  VW  antenna  at  Rome 
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Figure  29.  Optical  terminal  and  Rapidlink  antenna  at  Stockbridge 


4.6  Upgrade  -  Optical  System 

An  upgrade  to  the  optical  system  was  initiated  with  an  aim  towards  better  pointing  and  tracking 
performance  and  increased  communication  bandwidth.  Improvement  to  the  pointing  and  tracking  was 
needed  since  often  times  the  link  would  be  broken  by  external  factors  like  strong  winds  and  there  was 
no  automatic  mechanism  to  re-establish  the  link.  Manual  intervention  was  required  to  close  the  link 
again.  This  made  collection  of  long  term  link  performance  statistics  nearly  impossible.  A  new  control 
system  with  a  better  gimbal  promised  a  more  robust  pointing  and  tracking  system  which  would  enable 
better  data  collection  and  make  the  link  usable  for  quantum  communication  type  experiments. 

The  improvements  were  in  two  main  areas:  a  better  gimbal  from  Cobham  Sensor  Systems,  shown  below 
in  figure  29,  and  incorporation  of  its  coarse  pointing  functionality  into  the  system  control  loop.  The 
AOptix  terminal  takes  the  place  of  the  mounting  plates  shown  in  figure  29,  and  sits  between  the  gimbal 
forks.  The  general  idea  is  that  the  deformable  mirror  would  handle  the  correction  for  small  changes  in 
the  angle  of  arrival  of  the  signal  beam,  along  with  the  higher  order  corrections,  but  when  the  pointing 
moved  beyond  the  range  of  the  mirror,  the  control  system  would  tell  the  gimbal  to  make  a  gross 
correction  to  the  pointing.  Details  of  the  upgrade  follow.  In  the  fall  of  2013  AOptix  delivered  an 
upgraded  optical  system. 
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o 


Figure  30.  SPS-500  Positioner 


Laser  Terminals  Upgrade: 

1)  Control  circuit  electronics  and  software  upgrade  to  the  AOptix  laser  comm,  terminals  to  improve  link 
lock  and  performance. 

Specifically: 

PCBA:  PC  Board  Assembly 

AOC:  Adaptive  Optics  Control 

HVA:  High  Voltage  Amp 

WFS:  Wavefront  Sensor 

MM:  Membrane  Mirror 

NFOV:  Narrow  Field  of  View 

Upgrade  of  the  PCBAs  (AOC,  HVA,  WFS,  MM  electronics).  The  re-designed  AOC  PCBA  with  DSP1  and 
DSP2  upgrades  from  600  MHz  to  1  GHz,  improved  system  closed-loop  bandwidth  up  to  50  percent  (from 
l.OKhz  to  1.5Khz).  The  re-designed  MM  PCBA  facilitates  monitoring  of  environmental  conditions  of 
membrane  mirror  (MM).  Additionally,  the  new  circuit  design  supports  DWDM  reference  sources  at  1545 
and  1555nm  with  power  targets  of  0  to  20dBm.  The  re-designed  WFS  PCBA  circuit  improves  bandwidth 
matching,  noise  gain  peaking,  increased  dynamic  range  and  improved  sensitivity.  This  also  includes 
upgraded  software  for  the  system  that  allows  interaction  with  the  new  board  set,  among  other 
improvements. 
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Figure  31.  Previous  AOptix  control  Updated  AOptix  control  interface 

2)  Gimbal:  The  previous  FUR  (or  formerly  Directed  Perception)  PTU-D300  gimbals  were  not  designed  for 
continuous  pointing  in  the  control  loop.  They  are  for  rough  pointing  only.  Additionally,  the  PTU-D300 
had  a  pointing  resolution  of  0.006429°,  which  the  Cobham  Sensor  Systems  SPS-500  Positioner  and 
Control  System  improved  upon  with  ±  0.002°  accuracy.  Complete  SPS-500  gimbal  specs  are  in  Table  9. 


Table  9.  SPS  Gimbal  Specifications. 


Description 

Value 

Positioner 

Dimensions  (without  payload) 

23.41  H  x  18.23  W  x  8.0  D  inches 

Weight,  positioner/riser 

35.2  lb 

Weight,  payload 

35  lb  (estimated);  70  lb  (maximum) 

Weight,  total 

105.2  lb  (maximum) 

Azimuth  torque 

10  ±  1  ft  lbs. 

Elevation  torque 

5  ±  1  ft.  lbs. 

Azimuth  travel 

±  180° 

Elevation  travel 

-10°  to  +80°  ±  1° 

Mechanical  stops 

Azimuth:  ±  275°;  elevation:  -  35°  and  +  215° 

Azimuth  speed;  acceleration 

180°/sec;  180°/sec2 

Elevation  speed;  acceleration 

180°/sec;  180°/sec2 

Repeatability,  az  and  el 

±  0.002° 
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The  software  interface  differences  are  highlighted  in  figure  31,  showing  the  before  and  after  user 
interface,  respectively.  The  new  software  allows  user  control  of  the  adaptive  optics  implementation. 
Now  the  user  can  choose  between  full,  higher  order  correction  or  tip-tilt  only  correction.  This  can  be 
useful  in  demonstrating  the  improved  signal  performance  with  the  higher  order  correction. 
Additionally,  the  method  for  setting  the  best  loop  phase  for  the  feedback  control  has  been  automated 
with  a  simple  single  click  button  which  calculates  the  phase  number  and  sets  it.  In  figure  31  items  A-  D 
highlight  the  main  control  features  and  are  described  as  follows:  A  -  a  representation  of  the  32 
segments  of  the  wavefront  sensor  which  is  color  coded  by  the  number  of  detection  counts,  this  also 
corresponds  to  the  deformable  mirror  segmentation;  B  -  the  toggle  control  for  turning  the  control  loop 
on  or  off;  C  -  the  mode  for  setting  the  mirror  rest  state;  D  -  the  slider  bars  to  adjust  the  control  loop 
feedback  parameters. 

4.7  Results  -  Optical  Link 

The  optical  link  was  finally  re-established  from  Rome  to  Stockbridge  in  May  of  2014.  The  performance 
was  similar  to  that  of  the  initial  system  in  terms  of  actual  link  performance.  The  pointing  and  tracking 
control  system  has  not  been  fully  tested  due  to  many  systems  integration  issues  as  well  as  external  test 
site  issues.  An  example  of  the  measured  link  data  is  shown  in  figure  32  below.  The  drop  out  of  the 
signal  at  about  7:40  is  due  to  user  adjustments.  This  can  be  compared  to  the  previous  link  performance 
over  3  days'  time  from  May  2011  as  shown  in  figure  33. 
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Figure  32-  Link  performance  on  June  4th  2014 
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Figure  33.  Link  performance  over  24-26  May  2010 
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5  Photon  Entanglement  Transport  Through  a  Hyperspectral  Filter 

(Dr.  David  H.  Hughes  and  Dr.  Reinhard  Erdmann) 


5.1  Introduction 

Our  in-house  work  asked  the  question,  can  a  multi-access  hyperspectral  telescope  maintain 
entanglement  between  spontaneously  generated  photon  pairs  in  four  wave  mixing  processes,  thus 
providing  shared  secret  key  material  to  seed  stream  ciphers  for  data  encryption?  If  not,  can  the 
telescope  be  modified  to  do  so,  or  should  another  quantum  key  distribution  protocol  be  implemented 
that  does  not  rely  on  entanglement? 

Funded  by  AFOSR,  our  current  investigations  are  motivated  in  a  contractual  effort  funded  by  AFRL  6.2  to 
integrate  two  technologies,  a  keyed  quantum  communication  secret  key  distribution  system,  and  a 
multi-access  laser  communications  system.  The  former  technology  utilizes  entangled  photon  pairs  to 
distribute  secret  keys  and  is  highly  quantum  mechanical.  The  latter  technology  utilizes  a  Lyot  filter 
based  hyperspectral  tree  for  wavelength  division  multiplexing,  in  both  transmit  and  receive,  for  high 
data  rate  classical  optical  communications.  It  is  designed  for  communications  between  a  hub  and 
several  spokes. 


5.2  Background 

The  two  technologies  to  be  integrated  in  the  6.2  contractual  effort  are  shown  in  figure  34.  NuCrypt  LLC 
developed  the  entangled  pair  source  and  detection  systems  with  supporting  funding  by  an  AFRL/AFOSR 
STTR  [1],  while  Optical  Physics  Company  developed  the  multi-access  hyperspectral  telescope  under  an 
AFRL/SMC  SBIR  [2], 


Multi-Access  Lasercom 


Entangled  Pair  Generation  and  Detection 


Figure  34.  A  multi-access  hyperspectral  telescope  (left)  is  to  be  integrated  with  an  entangled  photon  source  and 

receiver  (right). 

The  hyperspectral  filter  configuration  on  the  left  of  figure  34  supports  four  simultaneous  users.  It  will  be 
modified  to  construct  three  telescopes,  each  supporting  two  simultaneous  users.  A  core  functional 
element  of  the  hyperspectral  telescope  is  the  filter  tree  comprised  of  Lyot  stages  [3].  Each  Lyot  stage 
comprises  a  Lyot  filter  stack  followed  by  a  polarizing  beam  splitter  (PBS)  made  of  silicon.  This  is  shown 
in  figure  35.  Incident  frequencies  chosen  from  the  ITU  grid,  separated  by  approximately  200  GHz,  are 
either  congruent  or  incongruent  with  respect  to  the  birefringence  of  each  Lyot  stack.  Basically,  the  Lyot 
stack  will  leave  invariant  the  polarization  state  of  a  congruent  frequency,  but  rotate  an  incongruent 
frequency  polarization  by  n/2.  The  polarization  beam  splitter,  rotated  by  tc/4  with  respect  to  the  Lyot 
stack  optical  axis,  projects  horizontal  and  vertical  inputs  with  respect  to  its  splitter  interface  into 
orthogonal  directions.  This  is  accomplished  by  total  internal  reflection  by  S  polarized  states,  here  the 
horizontal  states,  and  frustrated  internal  reflection,  or  transmission,  by  P  polarized  states,  here  the 
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vertical  states.  S  polarized  states  possess  polarization  states  orthogonal  to  the  plane  of  incidence.  P 
polarized  states  possess  polarization  states  parallel  to  that  plane;  they  tunnel  across  an  air  gap  at  the 
interface  between  the  two  silicon  prisms  comprising  the  PBS.  For  proper  operation,  all  frequencies  are 
chosen  from  the  ITU  200  GHz  grid.  We  denote  the  congruent  frequency  f1;  and  incongruent  frequency 

f2- 


Figure  35  Lyot  filter  stage  comprised  of  a  Lyot  filter  stack  followed  by  a  polarizing  beam  splitter,  or  PBS.  Here,  fi  is 
congruent  and  f2  incongruent  with  the  Lyot  stack.  Both  are  incident  at  7t/4  on  the  Lyot  stack,  but  f2  was  rotated  by 

7i/2  and  is  parallel  to  the  plane  of  incidence  of  the  PBS. 


a 

t 


Figure  36.  Entangled  photon  pair  amplitudes  enter  one  port  of  a  beam  splitter  (BS).  There  they  split  into  two 
paths.  They  recombine  when  entering  a  modified  Lyot  filter  stage  at  Port  1  and  Port  2;  classically,  we  would  say 

Port  1  or  Port  2. 

In  principal,  NuCrypt  LLC's  entangled  photon  pair  source  can  emit  entangled  pairs  collinearly  such  that 
those  pairs  arrive  along  a  line  at  the  Lyot  filter  stack.  That  will  not  be  implemented  in  the  6.2  integration 
project.  In  order  to  have  a  ghost  of  a  chance  at  accomplishing  a  quantum  key  distribution  protocol, 
classical  communications  are  necessary.  Leakage  between  the  classical  channels,  though  not  fatal  to  the 
classical  communications,  would  swamp  the  quantum  channel.  Instead,  in  the  6.2  contractual  effort  the 
quantum  and  classical  channels  are  split  right  away  and  processed  separately. 

Our  6.1  basic  research  thrust  here,  funded  by  AFOSR  LRIR13RI08COR,  focuses  on  modeling  photon 
entanglement  transport  through  a  modified  Lyot  hyperspectral  stage  in  the  absence  of  classical  signals. 
Shown  in  figure  36  is  the  modified  hyperspectral  stage.  Lyot  stacks  LI  and  L4  are  identical.  They  effect 
differential  changes  in  the  polarization  states  of  f2  and  f2  photons.  Lyot  stacks  L2  and  L3  are  identical  as 
well,  but  effect  no  polarization  changes  in  either  input  frequency.  Future  research  and  reporting  will 
vary  these  parameters.  For  now,  we  will  first  concentrate  on  just  LI  to  understand  the  behavior  of  the 
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entangled  state  transiting  it,  then  add  L4  and  a  beam  splitter  to  present  the  possibility  that  the  input 
entangled  state  enters  the  PBS  from  P2  or  P4. 


5.3  Entanglement  Transport 

Consider  an  entangled  photon  pair  input  into  the  Pi  port  of  the  Lyot  filter  stage  wherein  one  frequency 
is  congruent  with  the  Lyot  stack  and  the  other  is  incongruent.  Label  them  fx  and  f2  respectively.  They 
are  generated  in  a  spontaneous  four  wave  mixing  process  wherein  energy  and  photon  spin  are 
conserved.  A  quantum  joint  state  generated  by  the  process  can  be  written  in  terms  of  the  Lyot  stack 
reference  frame  or  the  polarization  beam  splitter  frame  (PBS)  as  below: 


la'">=V2 


(|f„h,P1),+|f„v,P,)1)  (|f2,h,P,)2+|f2,v,P,)2) 

- a - ® - Ti - 

(|f„h,P1),-|f„v.P,>1)  (|f2,h,P,)2-|f>,P,)2) 


^/2 


■J2 


(5.1a) 


=i(lf1.H.p1>1|L.H,P,)2+|f„v,pl)1|f2,v,pl>2) 


Symbols  fi  and  f2  refer  to  the  input  frequencies.  H  and  V  refer  to  horizontal  and  vertical  polarizations  in 
the  PBS  reference  frame,  while  symbols  h  and  v  refer  to  horizontal  and  vertical  polarization  with  respect 
to  the  Lyot  stack,  oriented  at  tc/4  with  respect  to  the  PBS  in  figure  35.  P2  is  the  leftmost  port  where  a 
state  first  enters  the  Lyot  stage,  as  shown  in  figure  35.  Commas  between  the  elements  within  a  ket  are 
shorthand  for  tensor  product  symbols,  product  signs  ensconced  within  circles.  A  comma  is  used  to  save 
space  on  the  page  within  an  equation.  Subscripts  on  the  kets  refer  to  the  fact  that  there  are  two 
photons  emitted  in  the  four  wave  mixing  process,  each  having  degrees  of  freedom  listed  within  the  ket. 
That  is,  each  photon  has  probability  amplitudes  associated  to  them.  Which  set  of  photon  amplitudes  to 
label  T  and  ‘2!  is  our  choice.  That  two  photons  are  emitted  by  a  physical  process  under  conditions  we 
engineer  a  possibility  to  achieve,  estimated  by  us  with  probabilities  of  occurrence,  is  not  up  to  us.  It  is 
decided  by  nature.  A  mathematical  way  to  express  this  is  to  say  that  each  photon  state  resides  in  their 
respective  Hilbert  spaces.  Those  Hilbert  spaces  are  entangled  or  they  are  not. 

We  will  often  write  (5.1a)  in  the  more  condensed  notation  of  (5.1b), 


=  ^(|f1,H,P1;f2,H,P1)  +  |fi,v,P1;f2,v,pi» 


(5.1b) 


Semicolons  divide  the  two  Hilbert  spaces.  Because  we  have  two  Hilbert  spaces  comprising  the  joint 
state,  the  subscripts  on  the  kets  sometimes  remain.  Position  wise,  however,  subspace  1  is  the  left  ket 
subspace  2  is  the  right  ket.  If  we  label  the  frequencies  pandv  those  subscripts  or  the  position  of  the 
kets  will  remain  1  and  2.  In  this  system,  frequency  is  immutable.  It  really  defines  the  existence  of  each 
Hilbert  subspace  in  the  joint  distribution.  In  essence,  it  is  not  so  hidden  a  variable. 

The  Lyot  stack  is  represented  by  a  series  of  unitary  operations  [6], 
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(3)  LiNbQ3  Plate 
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(5.2) 


In  the  Lyot  stack  frame,  the  Lyot  filter  is  the  product  of  the  unitary  operators  in  equation  5.2,  given  by. 
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Index  takes  values  {1,2},  and  refers  to  frequencies  fi  and  f2  .  <|).  is  the  optical  path  length  difference 
(OPD)  phase  for  each  plate  and  cp  is  the  cumulative  OPD  phase  for  the  stack.  Congruent  frequencies 
possess  cumulative  OPD  phases  cp  =  n,  for  a  total  OPD  phase  of  2n,  while  incongruent  frequencies 
possess  cumulative  OPD  phases  cp  =  n/2,  for  a  total  OPD  phase  of  n.  Thus  the  actions  of  the  Lyot  stack 
on  congruent  and  incongruent  frequency  states  is  a  function  of  the  input  frequency. 
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(5.4) 


where 


<*i  = 


'0  f 

'0  -i" 

'1  O' 

'1  O' 

5  ^9 

9 

,  and  crft  = 

1  0_ 

7  z 

i  0 

7  J 

.0 

7  u 

0  1 

(5.5) 


It  is  important  to  remind  oneself  that  the  Lyot  filter  stack  as  constructed  is  a  device  employing 
birefringence  to  operate  on  input  states  whose  frequencies  are  either  congruent  or  incongruent  with 
the  thickness  of  the  birefringent  plates  comprising  the  filter  stack.  This  allows  one  to  vary  the  input 
frequencies,  or  vary  the  birefringence  of  the  Lyot  filters  in  an  active  manner  to  gain  different  outputs,  a 
capability  about  which  more  will  be  reported  after  further  investigation. 

States  expressed  in  the  computational  basis  are: 


l»>- 


(5.6) 


So  that. 


,H|i>}  |v>] 
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,V2. 


4~2 


(|h)+|v))  and  |av)=[|h)  |v>] 
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jr(lhHV))  (5.7) 
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For  the  fixed  Lyot  filter  stacks  of  equation  (5.4),  their  action  on  the  congruent  and  incongruent  input 
states  of  equation  (5.1)  are,  suppressing  the  port  index  for  now. 


Lfi|ahi}  =  [|h),  Wjia, 

Lfi|avl)  =  [|h>1  |  v),  ]  i  cfj 


Lf2|ah2)  =  [|h)2  | v)2]ia2 

Lf2  Iav2>  =  [Ih>2  Iv)2]i°2 
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(5.8a-5.8d) 


But  for  a  phase,  the  Lyot  stack  leaves  invariant  polarizations  of  congruent  frequency  states,  while 
rotating  polarizations  of  incongruent  frequency  states  by  n/2.  The  two  frequencies  live  in  distinct 
Hilbert  spaces,  which  are  subspaces  of  the  joint  state.  That  fact  implies  unitary  operations  on  states  in 
one  subspace  do  not  effect  changes  in  states  in  the  other  subspace.  Both  states  propagate  collinearly 
through  the  Lyot  stack,  for  now  assumed  a  perfect,  lossless  device. 

Polarization  entanglement  in  the  joint  state  is  preserved  in  transit  through  the  Lyot  stack.  Prior  to 
entering  the  stack,  the  joint  state  given  in  equation  (5.1)  can  be  written  as  a  combination  of  two  Hilbert 
subspaces  via  a  Jordan-Schwinger  representation  [7,  8].  See  Appendix  5.A. 
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Or,  inserting  the  port  subscript  and  expressing  in  the  PBS  basis,  (5.1)  is  recovered. 


(5.9) 


[ l,  } =j=  (| : f,  >H.P,  >,  I  f2  •  H,  P,  }2  + 1  f, ,  V,  P, ),  |  f2 ,  V,P,  )2 )  (5.10) 


Under  the  action  of  the  Lyot  stack,  however,  and  according  to  equations  (4.8), 
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Or,  inserting  the  port,  PI,  and  using  the  condensed  notation, 

|aout)  =  ^(|f1,H,P1;f2,V,P1}  +  |f1,V,P1;f2,H,P1))  =  |bi„ 
|a0Ut)  =  ^(|fi5h,P1;f2,h,P1)  +  |f1,v,P1;f2,v,P1))  =  |bin) 


(5.11) 


(5.12) 


Equations  (5.12)  express  the  exit  state  into  PI  in  figure  35  in  terms  of  the  PBS  and  Lyot  reference  frames 
respectively.  The  input  state  to  the  system,  equation  (5.1)  is  rotationally  invariant.  The  output  state 
from  the  Lyot  stack,  'bin/  into  the  PBS,  is  not. 

Symbol,  R,  in  equations  (5.9)  is  the  operator  connecting  the  two  Hilbert  subspaces.  It  doesn't  always 
describe  an  entangled  state.  It  can  describe  all  four  of  the  so-called  maximally  entangled  Bell  states  if  R 
is  one  of  the  normalized  basis  elements  in  equation  (5.5).  The  Bell  states  span  the  four  dimensional 
space  of  two-qubit  composite  systems.  For  example,  a  linear  combination  of  them  can  describe  mixed 
states  that  are  not  entangled.  See  Appendix  5.A. 


Figure  37  Action  of  the  PBS  on  the  input  states.  A  vacuum  state  into  the  PBS  is  indicated  by  the  dashed  line.  We 
fill  that  vacuum  with  the  possibility  that  the  entangled  state  enters  at  P4  as  well  as  PI  by  placing  a  BS  in  the  input 

path  as  shown  in  Figure  36. 

We  now  add  the  possibility  of  the  state  entering  P4  as  well  as  PI,  as  shown  in  figure  35.  Two  photons  of 
an  entangled  pair  enter  the  BS  from  the  same  port.  They  both  exit  the  BS,  entering  PI  port  of  figure  36. 
The  modified  Lyot  stage  is  comprised  of  four  Lyot  filter  stacks  followed  by  a  polarizing  beam  splitter,  or 
PBS.  Here,  f2  is  congruent  and  f2  incongruent  with  the  Lyot  stacks,  LI  and  L4.  However,  in  this  model 
both  are  congruent  with  L2  and  L3.  M  is  a  mirror,  BS  is  a  non-polarizing  beam  splitter,  and  q  is  a  phase 
adjustment.  It  rectifies  the  phase  due  to  path  length  difference  and  reflections  between  the  two 
possible  ways  the  entangled  state  enters  the  Lyot  stage.  After  passing  through  the  BS,  the  state  input 
into  the  Lyot  stage  can  be  expressed  as, 
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(5.13) 


\  (I ' A  (HjP, ;  H2Pi  I  ;  ■ V2P, ))  +  |  ■ A  (HA ;  H2P4 1  V,P4 ;  V2P4 ))) 
|(|T0(H1P1;H2P4|V1P1;V2P4))+  |A |  V,P4;  V2P,))) 


The  'P  functions  are  called  Bell  state  forms.  They  collect  probability  amplitudes  possessing  different 
degrees  of  freedom  into  convenient  algebraic  forms.  See  Appendix  5.A.  Passing  through  the  Lyot  stack, 
and  then  the  PBS,  and  expressing  the  resulting  state  in  the  PBS  frame,  we  arrive  at. 


kut )  =  |  (I  (HA ;  V2P3  I  V,p3 ;  H2P2 ))  +  |  VP0  (HjPj ;  V2P2 1  V,P2 ;  H2P3 ))) 

|(|A(hA;v2p2|va;h2p3))  +  |^o(hip3;v2p3|v1p2;h2p2))) 


(5.14) 


Basically,  the  Lyot  stack  performs  rotations  on  the  input  state,  setting  it  up  for  path  projection  by  the 
PBS.  PBS  is  a  polarization  projector.  But  it  is  also  a  unitary  operator  [8,9].  (5.14)  contains  eight  states 
over  sixteen  single  constituent  bases.  Those  states  are  distributed  according  to  the  Bell  state  form  in 
which  they  are  contained.  Inspect  closely,  however,  the  degrees  of  freedom  in  both  (5.13)  and  (5.14.) 
As  sets,  they  are  all  different.  For  that  reason  every  one  of  the  above  forms  is  orthogonal  to  all  the 
others.  Since  each  Bell  form  is  normalized  by  the  square  root  of  two,  they  form  an  orthonormal  basis. 


Since  the  BS  breaks  the  symmetry  of  one  dimensional  collinear  propagation,  new  modes  appear.  Some 
have  very  desirable  properties.  Others  do  not,  at  least  for  our  purposes.  Here,  the  possibility  of 
simultaneous  arrivals  of  both  photons  to  one  or  the  other  party  occurs  on  average  half  the  time.  Now, 
that  might  be  useful.  It's  sharable  information,  but  it  conflicts  with  the  efficacy  of  another  protocol  that 
allows  unambiguous  information  sharing  between  the  legitimate  users.  To  ensure  unambiguous 
information  sharing  in  a  perfect  world  between  the  two  legitimate  users,  Alice  and  Ralph,  those 
measurements  are  discarded.  This  reduces  the  throughput  by  50%.  A  similar  situation  occurs  in 
traditional  BB84  and  Ekert91  under  basis  reconciliation. 


Appendix  5.B  outlines  a  transformation  on  5.14,  where  approximate  states  are  derived.  These  states  are 
those  arising  in  the  first  two  terms  in  (5.13)  and  then  (5.14).  The  transformation  of  the  state  output  by 
the  PBS,  truncated  and  renormalized  is  shown  below.  To  denote  that  they  arise  in  a  transformation,  we 
call  the  erstwhile  port/path  2,  path  'a/  and  port/path  3,  'b,'  and  constrain  the  angles  the  users  can 
measure  polarization  with  to  be  exactly  the  same.  The  common  angle  8  =  0a  +  0b  and  0a  =  0b.  The 
reduced  state  is, 


|dout)  =  -^^^(|A(ha;vbl  vb;ha))+|A(hb;va|va;hb))) 

+  sin(25)  /|^  (ha;hb  |  va;vb)^  +  |'F3  (hb;ha  |  vb;va))) 
v  2 

|dout)  =  -^^^-(1  ^SI95)  (1 1;00 1  00;1 1))+ 1 'fSI50)  (10;01 1 01;10>) ) 

+  S1^5)(|vy(228)  (i  1;10 1  01; 00)}  +  |'P(177)  (10;11 1 00;0l)} ) 


(5.15a) 


(5.15b) 
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(5.15c) 


|dj=  -^|^(|>Pl|,M)}+|>P|,i,)))  +  (|'Pf!>)  +  |'jf7> 


Defined  in  Appendix  5.A,  the  Bell  forms  are  not  unique  in  the  sense  that  they  were  constructed 
arbitrarily  by  collecting  terms  resulting  from  the  transformation  producing  (5.15).  The  amplitudes 
arising  in  that  transformation  do  split  up  into  distinct  Bell  form  classes,  but  the  amplitudes  within  a  class 
were  collected  and  summed  by  arbitrary  choice.  Be  that  as  it  may,  once  chosen  they  can  be  represented 
uniquely  by  binary  bit  strings,  which  can  then  be  expressed  in  base  ten.  Here,  h->l,  v->0,  a->l,  b->0. 
The  '|'  symbol  between  the  degrees  of  freedom  denote  +  or  -,  depending  on  the  particular  Bell  form  in 
which  they  are  ensconced.  In  this  condensed  notation, divides  the  two  subspaces  that  are  entangled. 
The  superscript  on  the  ¥  Bell  forms  is  the  base  10  representation  of  the  concatenated  degrees  of 
freedom  expressed  as  a  binary  stream.  Finally,  the  subscript  identifies  a  particular  Bell  form. 

An  interpretation  of  (5.15)  is  that  the  composite  system  is  behaving  like  a  single  system  possessing  four 
possible  states,  two  per  Bell  form  class,  and  all  transporting  entanglement.  In  the  absence  of 
decohering  processes,  this  obtains  even  when  the  two  spatial  modes  of  the  system  splay  the  parts  of  the 
system  extensively  over  space-time. 

To  compute  the  density  matrix  of  the  output  state,  we  encode  the  constituent  basis  elements  in  (5.15) 
by  their  internal  degrees  of  freedom. 
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Symbol  'c'  stands  for  cosine,  and  symbol  's'  for  sine.  This  8X8  matrix  contains  64  elements.  It  is  a 
reduced  dimension  version  of  the  full  system  matrix,  a  16X16  that  included  the  discarded  terms.  The  full 
system  matrix,  including  the  discarded  terms  is  displayed  in  Appendix  5.D.  Possessing  a  trace  equal  to 
one  (5.16)  is  a  proper  density  operator.  Moreover,  its  square  equals  itself,  so  it  describes  a  pure  state. 
Its  eigenvalues  are  the  set  {0,0,0, 0,0, 0,0,1} . 

One  measure  of  whether  a  state  is  entangled  is  obtained  by  forming  the  partial  transpose  of  a  density 
operator  purporting  to  represent  a  separable  state.  If  the  density  operator  of  a  two  particle  composite 
system  is  a  tensor  product  of  independent  density  operators,  and  each  density  operator  necessarily  has 
non-negative  eigenvalues,  as  do  their  transposes,  the  tensor  product  of  the  two  density  operators  will 
also  have  non-negative  eigenvalues.  Necessarily ,  so  will  the  density  operator  of  the  tensor  product  of 
one  system  with  the  transpose  of  the  other.  That  means  the  density  operator  is  separable  as  a  tensor 
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product  of  two  proper  density  operators.  However,  if  the  eigenvalues  of  the  partially  transposed 
density  operator  of  a  composite  system  has  at  least  one  negative  eigenvalue,  it  is  not  separable.  Its 
density  operator  cannot  be  represented  by  a  sum  of  tensor  products  of  two  qubits.  It  is  an  entangled 
system. 
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transpose  of 

(5.16) 

is  displayed  in 

(5.17). 

Its 

eigenvalues 

(5.17) 


The  partial 

{0,  0,0,0,  —1  /  2, 1  /  2, 1  /  2, 1  /  2}  .  Therefore,  because  at  least  one  eigenvalue  is  negative,  the  reduced 

constituent  amplitude  density  matrix  representation  of  the  composite  system  exiting  the  PBS,  (5.16), 
represents  an  entangled  system.  This  is  independent  of  angle  8. 


The  partial  transpose  technique  for  determining  separability  originated  with  Asher  Peres  [10],  and  its 
limitations  were  quantified  by  the  Horodeckis  [11,12].  Unfortunately,  the  method  is  not  scalable  in  the 
sense  that  some  density  matrices  of  higher  dimension  than  2X2  and  2X3,  known  to  be  non-separable,  do 
not  possess  any  negative  eigenvalues.  That  situation  has  been  termed  bound  entanglement.  The 
Horodeckis  determined  that  the  Peres  criterion  is  necessary  to  prove  non-separability,  but  is  not 
sufficient. 


It  is  interesting  to  trace  out  degrees  of  freedom  from  both  sub-systems  in  an  effort  to  ascertain  what 
degrees  of  freedom,  if  any,  are  left  entangled.  In  this  model,  the  frequency  degree  of  freedom  is 
associated  with  a  particular  ket  in  the  ordering  of  the  entangled  state,  as  discussed  below  (5.1).  Hence 
we  will  suppress  its  symbol  within  the  kets,  keeping  their  order  or  simply  denoting  the  internal  degrees  of 
freedom  associated  with  a  particular  frequency,  fl  or  f2,  by  '1'  and  '2'  subscripts.  Recalling  (5.15), 
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Its  composite  state  density  operator  can  be  written  as. 


cos 


Py  =■ 


|  \j/(195)  ^  ^\J/(195)  |  _|_  |  \J/(195)  ^  ^\J/(150)  |  _|_  |  \J/(150)  ^  ^\J/(195) 


+ 


xj/(15°)^vj/(150)  j 


+ 


cos(25)sin(25) 


f  |vf'|95)('F228|  +|'F|50 

+  |vI/228\(vFl95|  +|VF228 


'Ff8  +  PFj95 


'f‘77|+|vf;50)('F377|  a 


t;50  +  vf‘77  )(vf;95  +  vf/377)(vt/l50 


(5.19) 


+  Sm2^(  |'F3228)('Pf|  +|'F3228)('F"7|  +  |'F‘77}('F3228|  +  |^77)('F 


We  then  compute  the  partial  traces  over  the  different  degrees  of  freedom  within  the  states,  i.e., 
Tr^  x.  y)(p>j.)  and  Tr^x  y  The  first  is  a  partial  trace  over  path  and  the  second  a  partial  trace 

over  polarization. 

Tr(  x;  y)  (pT)  =  (aja2  |pT  |a1a2)  +  (a1b2  |pT  |a1b2)  +  (b1a2  |pT  |b[a2)  +  (b3b2 1 p^  jb^}  (5.20a) 

Tr(x_;y_)  (P*  )  =  (hlh2  |  Py  |  hlh2)  +  (hlV2  |  Py  |  V2)  +  (Vlh2  |  Py  |  Vlh2)  +  (V1V2  |  Py  |  V1V2  )  (5‘20b) 


Tracing  over  path  degrees  of  freedom  renders  the  reduced  density  matrix  in  polarization  degrees  of 
freedom,  (5.20a).  It  can  be  constructed  by  reference  to  (5.16)  and  look  for  terms  consistent  with  (5.20a). 
In  the  same  basis  used  in  (5.20a),  the  reduced  density  matrix  for  polarization  is. 
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We  have  written  the  binary  reductions  beneath  the  former  binary  bases,  having  traced  out  the  red  bits 
which  denoted  the  path,  leaving  polarization  degrees  of  freedom  intact.  Note  there  exists  just  four 
unique  bases  in  the  reduced  basis  set.  Two  rows  are  multiples  of  two  other  rows  in  (5.21).  This  causes 
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degeneracy  in  the  eigenvalue  spectrum.  Eigenvalues  of  (5.21)  are  the  set  {0,0, 0,0, 0,0, 1/2, 1/2}.  Its  partial 
transpose  is 
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Matrix  (5.22)  has  two  negative  eigenvalues  in  the  set  {-1/4,  -1/4,  1/4,  1/4,  1/4,  1/4,  1/4,  1/4}.  Hence  the 
reduced  density  in  polarization  necessarily  describes  an  entangled  state.  It  reduces  to  (5.23),  on  the  left. 
Its  partial  transpose  is  on  the  right.  Respectively,  their  eigenvalues  are  {0,0, 0,1}  and  (-l/2,l/2,l/2,l/2}. 
Compare  these  eigenvalues  with  those  found  for  (5.16),  (5.17),  and  for  the  full  system  density  matrix  in 
Appendix  5.C. 
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prp0i  is  pure,  for  it  equals  its  square.  The  partial  transpose  has  one  negative  eigenvalue.  Polarization 
entanglement  is  maintained,  and  over  all  common  measurement  angles  performed  by  nodes  a  and  b. 

Returning  to  (5.16)  we  now  trace  over  the  polarization  degrees  of  freedom  according  to  the  prescription 
in  (5.20b). 
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partial  transpose  of  (5.24)  is. 
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(5.24) 


(5.25) 


In  the  current  configuration,  the  reduced  density  matrix  for  path  equals  its  partial  transpose.  Their 

moon  cos2 (25)  cos2 (25)  sin2 (25)  sin2  (28), 

eigenvalues  are  {0,0, 0,0, - , - , - , - j.  These  eigenvalues  are  always 


positive.  If  we  put  the  frequency  degrees  of  freedom  explicitly  back  into  the  bases,  though  it  was  always 
there  implicitly  in  the  ordering  and  identification  of  the  kets  as  1  and  2,  nothing  changes. 
Frequency/path  is  still  a  mixed  state  in  its  reduced  density  after  the  polarization  has  been  removed. 
Correlations  exist,  but  not  entanglement  between  those  degrees  of  freedom.  The  Lyot  stack  acts  on 
frequency  states  to  either  change  their  polarization  or  leave  it  effectively  invariant.  The  PBS  projects 
those  states  into  different  paths,  depending  on  their  polarization.  Under  the  current  configuration, 
polarization  is  the  entanglement  carrier  within  the  extant  degrees  of  freedom.  That  does  not  necessarily 
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mean  that  path  and  polarization  and  frequency  are  not  entangled  in  the  joint  distribution.  (5.16)  proves 
otherwise. 


Returning  to  (5.25),  we  see  that  there  exist  just  two  unique  basis  elements:  (a;b)  and  (b;a).  This  means 
that  frequency  fx  will  appear  at  port  (a)  P2,  or  (b)  P3,  with  equal  probability.  Similarly,  frequency  f2 
appears  at  either  port.  But  that  coincidence  does  not  mean  that  path  and  frequency  are  entangled.  They 
are  correlated,  but  not  entangled  in  the  reduced  matrix  representation  of  the  system.  Moreover,  the 
reduced  density  matrix  for  path  and  frequency  does  not  equal  its  square.  It  is  not  a  pure  state.  However, 
had  we  included  the  full  set  of  possibilities,  the  two  photon  coincidence  states,  path  and  frequency 
remain  entangled,  i.e.. 
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(5.26) 


(5.26)  is  a  pure  state,  and  its  partial  transpose  does  possess  a  negative  eigenvalue. 


Another  measure  of  entanglement,  different  from  the  negativity  in  the  partial  transpose  eigenvalue 
spectrum  of  the  constituent  particle  system  density  matrix,  is  the  comparison  of  probabilities  one  would 
get  if  the  system  were  classical  and  separable  with  those  probabilities  computed  using  quantum 
mechanics.  The  concept  was  introduced  by  John  Bell  [13]  in  response  to  Einstein,  Podolsky,  and  Rosen's 
conjecture  that  quantum  mechanics  is  incomplete;  it  did  not  explain  correlations  between  non¬ 
interacting  particles  possessing  observables  simultaneously  measured  in  an  inertial  reference  frame  and 
the  two  parties  separated  by  space  like  metric  distances  [14].  The  particles  in  question,  EPR  argue, 
possess  real  properties  even  in  the  absence  of  measurement,  and  that  no  connection  between  the 
particles  can  exist  that  can  propagate  at  speeds  greater  than  that  of  light  in  vacuum  such  that 
measurements  on  one  affect  outcomes  of  the  other.  Bell's  theorem  argues  against  the  validity  of  EPR's 
assumptions  on  local  realism,  that  quantum  systems  possess  definite  values,  and  that  quantum 
mechanics  predicts  the  outcome  of  correlation  experiments,  whereas  classical  statistical  mechanics 
cannot.  His  concept  has  been  extended  by  several  researchers,  such  as  Clauser,  Horne,  Shimony,  and 
Holt  [15]  who  formulated  the  CHSH  inequality.  See  Appendix  5.C,  where  we  compute  the  so-called  CHSH 
inequality  following  a  prescription  given  by  Stenholm  and  Suominen  [16]. 


5.4  Quantum  Key  Distribution 

Returning  now  to  (5.15a),  we  expand  it  to  make  explicit  the  degrees  of  freedom  and  their  structural 
relationship. 
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=  -^Pd^,  (1 1;00 1  00;ll))+|'P1  (10;01 1  01;10)) ) 
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(5.27b) 


(5.27c) 


Choosing  the  PBS  basis,  8=0,  and  having  agreed  to  discard  two  photon  coincidence  per  node  events, 
polarization  measurements  at  'a'  and  'b'  are  100%  anti-correlated.  Choosing  the  alternative  basis,  8=p/4, 
polarization  measurements  are  100%  correlated. 


Polarization  measurements  performed  by  parties  at  'a'  and  at  'b',  call  them  Alice  and  Ralph,  yielding  two 
photon  events  can  be  discarded.  That  amounts  to  an  average  50%  reduction  in  throughput.  It  happens 
in  BB84  when  Alice  randomly  chooses  a  basis  from  two,  two  dimensional  bases  in  which  to  send  a 
photon,  and  Ralph  independently  and  randomly  measures  in  one  of  the  two,  bases.  They  retain  only 
those  bases  that  coincide,  on  average  50%  of  the  time.  Here,  however,  even  when  they  have  chosen  to 
stick  with  two  possible  bases  differing  by  45°,  they  would  still  have  to  reconcile  their  choices.  That  is,  if 
they  choose  those  bases  at  random  and  independently,  half  the  time  their  choices  would  agree  and  half 
the  time  not.  That's  a  total  reduction  in  throughput  to  25%.  This  can  be  avoided.  If  they  share  a  secret 
key  allowing  them  to  select  predetermined  bases,  basis  reconciliation  is  unnecessary  [17].  That  key  can 
be  generated  by  a  classical  stream  cipher  possessing  enormous  mathematical  complexity.  To  generate 
long  running  keys,  all  it  needs  is  a  secret  seed  key.  That  might  be  obtainable,  at  least  in  principle,  by  QKD. 


Table  10  summarizes  the  ideal  results  two  possibly  displaced  users  would  measure  on  the  composite 
state  after  discarding  two  photon  arrivals  per  port,  and  when  the  input  state  utilized  is, 


|ain>  =  |(|H1P1;H2P1+V1P1;V2P1)+  |H,P4;H2P4 1  V.P,;' V2P4» 


(5.28) 


Table  10  Polarization  Measurements  by  Two  Displaced  Users 


5  =  0 

5  =  tt/4 

Alice 

H 

V 

h 

V 

Ralph 

V 

H 

h 

V 

5.5  Conclusion 

Designed  for  wavelength  division  multiplexing  of  classical  electromagnetic  optical  beams,  the  Lyot  stack 
employs  birefringence  to  discriminate  between  frequencies  of  incoming  states.  Their  optical  path  length 
(OPD)  phases  differ  when  passing  through  an  engineered  thickness.  Those  frequencies  whose  OPD 
phase  is  an  integral  multiple  of  2n  upon  exit  from  the  stack  are  termed  congruent  frequencies. 
Frequencies  whose  OPD  phase  is  an  odd  integral  of  n  are  termed  incongruent  frequencies.  Polarization 
states  of  congruent  frequencies  are  left  invariant,  while  those  of  incongruent  states  are  rotated  by  n/2. 
After  exiting  the  Lyot  stack,  the  two  frequencies  are  split  into  different  output  paths  for  further 
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processing.  This  allows  multiple  access  to  a  hub  whose  spokes  are  geographically  displaced.  An  obvious 
application  is  reception  and  dissemination  of  ISR  data  collected  by  the  spokes,  sent  to  the  hub,  and 
transported  either  by  satellite  cross  links  or  by  direct  linkage  from  satellite  to  another  location  not 
accessible  by  the  spokes. 

Our  in-house  6.1  basic  work  thus  far  focused  on  attempting  to  test  the  hypothesis  that  polarization 
entanglement  of  a  rotationally  invariant  entangled  quantum  state  is  not  preserved  in  transit  through  a 
perfect  birefringent  spectral  filter  stage.  In  the  current  model,  our  analysis  has  shown  that,  with 
modification,  polarization  entanglement  is  preserved  in  the  joint  state  as  it  transits  through  a  perfect 
birefringent  spectral  filter  stage  comprised  of  a  Lyot  stack  followed  by  a  polarization  beam  splitter. 

Initial  polarization  entanglement  is  maintained  through  the  Lyot  stack,  but  is  changed  from  a  rotationally 
invariant  entangled  state  to  a  non-rotationally  invariant  state.  Passage  through  the  PBS  still  maintains 
polarization  entanglement  in  the  joint  state.  In  the  joint  state,  polarization  is  entangled  with  path.  One 
modification  ensuring  this  is  a  purification  operation  allowing  the  possibility  that  the  state  enters  not 
one,  but  two  ports.  This  is  accomplished  physically  by  placing  a  non-polarizing  beam  splitter  in  the 
entrance  path  to  the  Lyot  stage.  The  state  exiting  the  PBS  is  then  transformed  to  a  linear  combination  of 
pure  states  allowing  users  to  execute  polarization  measurements  in  either  the  PBS  basis  or  the  Lyot  filter 
basis  rotated  by  45°.  This  capability  comes  at  a  cost.  Cross  terms  from  the  BS  transformation  also  enter 
the  Lyot  stage,  and  they  produce  the  possibility  of  two  photon  coincidence  events  at  each  PBS  output 
port.  If  the  users  discard  those  two  photon  coincident  events,  polarization  measurements  in  the  PBS 
basis  are  100%  anti-correlated,  while  measurements  in  the  Lyot  stack  basis  are  %100  correlated.  This 
incurs  on  average  a  50%  reduction  in  throughput.  Of  course,  the  users  may  wish  to  use  the  discarded 
events. 

Experimental  results  in  the  literature  indicate  path  and  frequency  can  become  entangled  [18].  If  the  BS 
is  removed  from  the  collinear  input,  only  one  entry  port  is  available.  After  passing  through  the  Lyot 
stack,  fl  and  f2  randomly  exit  either  P2/(a)  or  P3/(b),  but  always  at  different  ports.  They  do  so  because 
the  exit  state  from  the  PBS  is  |HP2;  VP3>+|VP3;  HP2>;  polarization  H  always  goes  out  P2,  and  polarization 
V  always  goes  out  P3.  Frequencies  associated  with  the  particular  polarization  states  follow.  A 
measurement  of  frequency  allows  Alice  to  infer  Ralph's  frequency.  On  that  basis,  information  is  shared. 
How  securely  that  information  is  shared  will  likely  require  a  buttressing  protocol,  perhaps  involving 
shared  time  bin  detection  [19]. 

A  temporal  delay  will  occur  when  an  eavesdropper  intercepts,  then  detects  the  frequency  for 
exploitation,  before  sending  that  photon  on  its  way  to  the  legitimate  user.  Accurate  and  precise  time 
binning  might  be  able  to  detect  that  temporal  delay  and  indicate  the  intrusion  compromising  the 
confidentiality  of  the  shared  information.  In  that  case,  one  needs  only  to  flip  the  logic  value  assigned  to 
the  frequency.  That  simple  operation  imparts  a  coin  flip  guess  on  its  value  to  the  eavesdropper. 
Furthermore,  there  is  no  need  to  maintain  a  common  rotation  angle,  8,  for  polarization  measurements. 
Availability  is  more  robust. 

Simpler  methods  exist  for  this  type  of  information  sharing,  protocols  preferred  by  engineers  tasked  to 
get  the  job  done  with  optimal  information  confidentiality,  integrity,  and  availability.  Our  interest  in  the 
current  Lyot  stage  device  characteristics  goes  beyond  QKD.  A  future  report  will  convey  details. 
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Appendix  5.A  Generalized  Bell  Forms 

A  general  bipartite  combination  of  two  states  in  the  Jordan-Schwinger  representation  is  given  by, 
|4')=[|x1)|w1)]R.a 

=  (R0+R3)|Xl)|y2)  +  (Rl-iR2)|Xi>|Z2)  +  (Rl+iR2)|wi>|y2>  +  (Ro-R3)|wi>|z2) 

^  V^/R0,R1,R2,R3  { {  X 1  ’ W 1 }  5  {  y 2  ’  Z2  } } 

For  R0  =  -j=  and  Ri  =  R2  =  Rb  =  0,  and  the  degrees  of  freedom  carried  by  the  state,  (Xj;z2  |  wty2 
where  the  semicolon  differentiates  between  the  two  spaces, 

=  ^(1 X1 ) I y2 >  + 1 Wi ) | z2 ))  □  |^0) 

Similarly, 

XF0,Jf,0,o)=^(lXl)lZ2>  +  |Wl)|y2))  D  W 
XFo,0,-j=,0^  =^(|Xl)|Z2>-|Wl)|y2))  D  |^2> 

XF°,°,0^)=  ^(lXl)ly2)"lWl)lZ2))  D  W 

In  the  Lyot  stack  reference  frame,  the  Lyot  stack  changes  the  entangled  (5.A.2)  state  to  the 
entangled  (5.A.5)  state.  The  transformation  in  Appendix  5.b  recovers  this  last  as  part  of  the 
pure  state  derived  there  as  a  linear  combination  of  Bell  state  forms  in  the  PBS  frame  and  the 
Lyot  frame. 


|y2> 

I  Z2  ) 


=  [K)  h)] 


ro+r3  Rj-iR2 

Rj+iR2  R0-R3 


|y2) 


L1R1  I 

-4, 0,0,0 

x/2 


^"o  /  *\  f  •  ^3 
~  1  r—  (“0^2  —  — 1  nr 


V2 


Hence,  in  the  Lyot  stack  frame,  diagonal  to  the  PBS  frame, 


(5.A.  1) 

(5.A.2) 

(5.A.3) 

(5.A.4) 

(5.A.5) 

(5.A.6) 
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'i'H|hl}|v,)]L1RL  L’ 


h2> 


=  [Ihl)  |Vl)] 


_,0,0,0 

72 

-i-7=  0 

V2 


0  i 


ll. 


K) 


=  -^(lhi»lh2>-|v1)|v2» 


(5.A.7) 


The  state  is  polarization  entangled. 

Traditional  Bell  states  usually  only  consider  one  degree  of  freedom.  The  Bell  forms  used  in  this  report  are 
orthogonal,  even  within  the  same  Bell  form  but  carrying  different  degrees  of  freedom.  Thus  these  Bell 
forms  are  generalizations  of  traditional  Bell  states. 
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Appendix  5.B  Transformation  of  the  Composite  State 


Figure  38  Two  coordinates  systems  related  by  a  tt/2  rotation  about  a  common  axis,  with  local  coordinate  rotations 
about  their  respective  propagation  directions.  P2->  'a'  and  P3  ->  'b'  in  the  transformed  frame. 

For  convenience,  we  will  use  very  compact  notation.  A  single  state  ket  and  a  bipartite  ket  will  be 
specified  in  the  following  ways. 


|'F1)  =  |f,h,a)1  =|h‘)  |vF12)  =  |fl,h,a)i|f2,v,b)2  =  |  f  1,  hf )  |  f  2,  )  =  |  hf ;  > 

In  the  PBS  frame,  with  P2  — >  a  P3  — >  b ,  the  state  exiting  the  PBS  is, 

I  cout )  =  j  (|  %  (f  1 ,  (Hja,  Vjb);  f  2,  (H2a,  V2b)))  +  |  T,  (f  1,  (Hjb,  V,a);  f  2,  (H2b,  V2a)))) 
+  j(|^(fl, (Hia, V,b); f  2, (H2b, V2a)))  +  |  ¥,  (f  1 , (Hjb, V.a);  f  2, (H2a, V2b)))) 


(5.B.1) 


(5.B.2) 


Note  that  product  of  the  frequencies  are  homogeneous  throughout,  in  fact  are  implicit  in  the  labeling  of 
the  kets  1  and  2.  Then  (5.B.2)  is  written. 


lco,)  =  |(|fl),®|f4)  (|>r,((Hla,Vlb);(H,a,V2b))>+|'P]((Hlb,Vla);(H,b,V2a))» 
+i(|fl>|®|f2)2)(|4'1((H,a,V1b);(H2b,V2a)))+  |^,((H,b,V1a);(H2a,V2b)))) 


(5.B.3) 


Performing  arbitrary  independent  rotations  at  both  output  ports  to  perform  measurements  of 
polarization  at  two  displaced  positions  upon  receive  gives. 
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Id)  =  Wf  (|hi;h2)  +  |hi;h2)_|Vi;V2)_|Vi;V2))  Sin(0a+0b) 

+  2 (  K’V2)+  Vl’;ll2)  +  |lll>;V2)  +  |Vnll2))  COS(0a  +  0b  ) 

+^f  (|K;v^)  +  |v®;h*))  cos(2eB)  +  ^(|hf;v$)  +  |v|,;h$))  cos(20b) 
+lif  (lhi;h2)"|vi’v2))  sin(20a)  +  ^f  (|hi;h2}-|vi;v^))sin(20b) 


(5.B.4) 


We  switched  to  lower  case  here  to  denote  we  are  working  in  the  transformed  frame.  Note  the 
occurrence  of  photons  simultaneously  arriving  at  either  'a'  or  'b'.  These  arise  from  the  mixed  terms 
exiting  the  BS.  The  consequence  of  this  two  photon  coincidence  effect  is  that  it  presents  noise  to  a 
particular  QKD  protocol  relying  on  100%  correlation  in  one  basis,  and  100%  anticorrelation  in  a  45° 
rotated  basis. 

Though  not  necessary,  If  0  =  0b  and0  +0b  =5,  (5.B.4)  reduces  to. 


|h^;h5)-|vr;v5)  +  |hIb;h;)-|vIb;^)^ 
275  yh-h^-|v-va2)  +  |h?;h^-|v^v^) 

_^nhi;V2)  +  |Vl’h2)  +  |hl’V2)  +  |Vi;h2) 

7 ""  v+  h; ;v2^+  Vj;h2^+ 


jin  (25) 


(5.B.5) 


Combining  the  degrees  of  freedom  within  the  kets  into  Bell  forms,  and  encoding  them  as  binary  strings. 


Y3(l  1;10 1  01; 00)  +  AF3(10;1 1 1  00; 01)  + 

¥3  (1 1;  1 1 1 0 1;  0 1)  +  ¥3  (1 0;  1 0 1 00;  00) 


1  f'F1(ll;00|00;ll)  +  'F1(10;01|01;10)+^| 

+  —  cos  25 

2  ^'F1(l  1;  0 1 1  0 1;  1 1)  +  'P1(10;00 1 00;10)  J  v  ’ 


(5.B.6) 


Then,  taking  the  entire  string  of  degrees  of  freedom  without  the  ‘\’  and  the  place  the  base  10 
equivalent  of  the  result  as  a  superscript  to  express  the  states  as. 


=  ^  (V,228)  +  vF(,I77)  +  'pf40  +  'P(3160))sin(25) 
+  ^  +  'P?50*  +  ^!215)  +  ^!‘30))  cos (25) 


(5.B.7) 


This  saves  space  and  uniquely  expresses  the  distribution  of  the  degrees  of  freedom  by  the  Bell  forms. 
The  last  two  terms  in  each  state  are  cross  terms  exiting  the  BS,  inducing  two  photon  coincidences  per 
port.  If  the  other  terms  are  projected  out,  meaning  the  measurements  on  the  cross  terms  are  removed 
from  consideration,  the  state  of  interest  for  the  current  application,  a  particular  QKD  protocol,  is  found  as 
follows.  Define  a  projector  onto  the  states  of  interest. 
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where  j  ={1,3}  and  d={228,177,195,150}.  Then  operate  on  |cout)  and  renormalize,  i.e., 

|p\_  A|d)  _  ^Vf/O95)  |  COS(28)  / w(228)  vt,(177)\  sin(28) 

i  j  >/2  1  3  3  j  ^ 


(5.B.9) 


(5.B.9)  means  Alice  and  Ralph  utilize  only  those  measurements  involving  the  states  entering  as  a  whole 
through  PI  and  P4  in  figure36  and  disregard  the  rest. 
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Appendix  5.C 


The  CHSH  inequality  states  that  the  maximum  classically  achievable  value  is  two.  Following  closely  an 
example  given  in  [15],  a  Bell  test,  or  BT,  is  the  CHSH  inequality  given  in  equation  (5.C.I.)  It  is  the 
relationship  between  expectation  values,  or  mean  values  of  measurements  made  by  'a'  and  'b'  on  some 
observable,  perhaps  the  phase  induced  on  polarization  states  by  a  phase  plate.  Classically, 

BT  =  |  E(a1,b1)  +  E(a2,b1)  +  E(a1,b2)-E(a2,b2)  |  <2  (5.C.1) 

Two  parties  perform  ensemble  measurements  on  two  observables  each.  Then  they  compute  the 
expectation  values  and  insert  their  results  into  (5.C.1).  The  maximum  classical  value  of  BT  is  2,  because 
the  normalized  expectation  values,  computed  by  a  classical  probability  density  function,  each  achieve  a 
value  of  1.  If  BT  is  violated  by  computing  the  expectation  values  using  quantum  mechanical  probability 
amplitudes,  quantum  correlations  exist  between  the  entangled  constituents  in  the  probability  amplitudes 
describing  the  state  of  the  system. 

Consider  observables,  Oa((|)a)  and  Ob  (<|)b),  defined  as  follows, 

0  ei<t>a  1  T  0  e^b~ 

0(6)=  and  Oh(d>„)=  (5.C.2) 

aVVa;  e-i<t>a  Q  e_4b  Q  V  ' 

-  Ja  L  Jb 

They  act  only  on  the  polarization  degrees  of  freedom  within  each  Bell  form.  Nodes  'a'  and  'b'  make 
measurements  on  (j)a  and  (j)b  respectively.  Suppose  'a'  and  'b'  have  agreed  to  perform  their 

measurements  in  their  respective  diagonal  bases,  i.e.,  where  the  observed  polarizations  are  100% 
correlated.  Using  the  computational  basis  to  express  the  polarization  states,  the  expectation  values  of 
interest  are. 


(4.C.3)  reduces  to  (4.C.4). 

E(<l>aA)=-COS(4>a  +  4>b)  (5-C-4> 


The  CSHS  inequality  is  violated.  The  state  possesses  non-local  correlations. 
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All  cross  class  Bell  form  terms  are  zero  in  (5.C.11),  leaving  only  the  terms  computed  above.  Those 
surviving  terms  are  weighted  by  the  sinusoidal  factors.  Then  the  expectation  value  for  Oab  is, 

E(<|>a,<|>b)=(doUt  |Oab  | dout )  =  cos2  (S)cos(^a  -  <|>b)  +  sin2  (8)cos(^)a  +  <|>b)  .  (5.C.11) 


=  2V2  (cos2  (28)  +  sin2  (28))=  2V2  (5.C.12) 

(5.C.12)  is  independent  of  angle  8 .  We  emphasize  that  cross  terms  from  the  BS  giving  rise  to  two 
photon  coincidences  in  either  path  have  been  discarded. 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 

58 


Appendix  5.D 


The  full  system  matrix  is  given  in  Table  11.  It's  the  matrix  on  the  left.  On  the  right  is  its  partial  transpose. 
Their  eigenvalues  are  the  sets  {15  zeroes,  1}  and  {12  zeroes,  -1/2, 1/2, 1/2, 1/2}  respectively. 
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6  Summary 

Title:  Investigating  Quantum  Data  Encrypted  Modulation  States 

Reporting  Period:  FY14 

Lab  Task  Manager:  Dr.  David  H.  Hughes 

Phone:  315-330-4122/DSN  587-4122/  Email:  david.hughes.16@us.af.mil 

AFOSR  Program  Manager:  Dr.  Arje  Nachman 

Research  Objectives:  Investigate  modulated  optical  signals,  both  quantum  and  classical,  for 
information  confidentiality,  integrity,  and  availability 

6.1  Technical  Work 

Simulation  and  measurement  of  the  noise  sources  in  coherent  states  of  light  is  being  conducted.  Over 
and  above  the  quantum  noise  of  light,  in  which  recent  physical  cryptographic  systems  utilize  to  impart 
randomness  to  a  ciphertext,  additional  noise  sources  can  have  an  even  more  debilitating  effect  on 
attempts  to  decrypt  ciphertext.  These  latter  sources  of  noise  include  detection  system  thermal  noise  as 
well  as  source  laser  phase  drift  noise. 

A  basic  comparative  study  between  long  wavelength  infrared  (LWIR)  verses  short  wavelength  infrared 
(SWIR)  propagation  in  fair  and  inclement  weather  is  being  conducted.  Currently  the  communication 
bands  under  study  are  1.5  and  10  micron  bands.  Despite  the  huge  infrastructure  devoted  to  the  1.5 
micron  band,  work  in  the  10  micron  band  by  several  workers,  including  AFRL/RITE  engineers,  indicates 
that  propagation  in  the  10  micron  band  is  more  robust  in  free  space  optics.  That  translates  to  greater 
information  availability. 

A  combined  optical/RF  network  is  under  development  between  Stockbridge  and  Rome.  Currently,  the 
optical  band  is  primarily  1.5  micron  band  carrier  frequencies;  turbulence  mitigation  is  accomplished  by 
curvature  adaptive  optics.  Planned  is  deployment  of  the  LWIR/SWIR  comparison  experiment  system 
over  the  longer  range  provided  by  the  Rome/Stockbridge  link  as  well  as  a  coherent  state  quantum  data 
encryption  system. 

We  asked  the  question:  Is  polarization  entanglement  preserved  in  transit  through  the  device?  If  not,  can 
the  device  be  modified  to  preserve  the  entanglement?  Through  modeling  we  found  indications  that 
polarization  entanglement  can  be  preserved  and  used  in  quantum  key  distribution. 


Table  12  Resource  Dispersal  FY14 


Government  Salaries  (6.1)  Basic 

$107,000 

Onsite  Contractor  Support  (6.1  Basic) 

$99,640 

Equipment  (6.1  Basic) 

$73,360 
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